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ABSTRACT 

We present source catalogs for the 4 Ms Chandra Deep Field-South (CDF-S), which is the deepest Chandra 
survey to date and covers an area of 464.5 arcmin^. We provide a main Chandra source catalog, which contains 
740 X-ray sources that are detected with WAVDETECT at a false-positive probability threshold of 10"^ in at least 
one of three X-ray bands (0.5-8 keV, full band; 0.5-2 keV, soft band; and 2-8 keV, hard band) and also satisfy 
a binomial-probability source-selection criterion of P < 0.004 (i.e., the probability of sources not being real 
is less than 0.004); this approach is designed to maximize the number of reliable sources detected. A total of 
300 main-catalog sources are new compared to the previous 2 Ms CDF-S main-catalog sources. We determine 
X-ray source positions using centroid and matched-filter techniques and obtain a median positional uncertainty 
of w 0.42". We also provide a supplementary catalog, which consists of 36 sources that are detected with 
WAVDETECT at a false-positive probabiHty threshold of 10"^, satisfy the condition of 0.004 < P < 0. 1, and have 
an optical counterpart with R < 24. Multiwavelength identifications, basic optical/infrared/radio photometry, 
and spectroscopic/photometric redshifts are provided for the X-ray sources in the main and supplementary 
catalogs. 716 (w 97%) of the 740 main-catalog sources have multiwavelength counterparts, with 673 (sa 94% 
of 716) having either spectroscopic or photometric redshifts. The 740 main-catalog sources span broad ranges 
of full-band flux and 0.5-8 keV luminosity; the 300 new main-catalog sources span similar ranges although 
they tend to be systematically lower Basic analyses of the X-ray and multiwavelength properties of the sources 
indicate that > 75% of the main-catalog sources are AGNs; of the 300 new main-catalog sources, about 35% 
are likely normal and starburst galaxies, reflecting the rise of normal and starburst galaxies at the very faint 
flux levels uniquely accessible to the 4 Ms CDF-S. Near the center of the 4 Ms CDF-S (i.e., within an off-axis 
angle of 3'), the observed AGN and galaxy source densities have reached 9800^} jqq deg"^ and 6900;!;goQ' deg"^, 
respectively. Simulations show that our main catalog is highly reliable and is reasonably complete. The mean 
backgrounds (corrected for vignetting and exposure-time variations) are 0.063 and 0.178 count Ms"' pixel"' 
(for a pixel size of 0.492") for the soft and hard bands, respectively; the majority of the pixels have zero 
background counts. The 4 Ms CDF-S reaches on-axis flux limits of « 3.2 x 10"''^, 9.1 x 10"'^, and 5.5 x 10"''' 
erg cm"^ s"' for the full, soft, and hard bands, respectively. An increase in the CDF-S exposure time by a factor 
of ss 2-2.5 would provide further significant gains and probe key unexplored discovery space. 

Subject headings: cosmology: observations — diffuse radiation — galaxies:active — surveys — X-rays: galax- 
ies 
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1. INTRODUCTION 

Deep X-ray surveys indicate that the cosmic X-ray back- 
ground (CXRB) is largely due to accretion onto supermassive 
black holes integrated over cosmic time. One of the greatest 
legacies of the Chandra X-ray Observatory {Chandra) is the 
characterization of the CXRB sources thanks to its extraordi- 
nary sensitivity. The Chandra Deep Field-North and Chan- 
dra Deep Field-South (CDF-N and CDF-S, jointly CDFs) 
are the two deepest Chandra surveys (see Brandt & Hasinger 
2005 and Brandt & Alexander 2010 for reviews of deep ex- 
tragalactic X-ray surveys), each covering « 450 arcmin^ areas 
with tremendous multiwavelength observational investments. 
Most of the CDF sources are active galactic nuclei (AGNs), 
often obscured, at z « 0.1-5.2. The CDFs have found the 
highest density of reliably identified AGNs on the sky, with 
an AGN source density approaching ten thousand sources per 
deg^ (e.g., Bauer et al. 2004). At faint fluxes, the CDFs are 
also detecting large numbers of starburst and normal galax- 
ies at z w 0.1-2 as well as a few individual off-nuclear X-ray 
binaries at z « 0.05-0.3. 

Deeper X-ray observations not only further improve the 
photon statistics that are required to understand better the al- 
ready detected sources via X-ray spectral and variability con- 
straints, but also probe further down the X-ray luminosity ver- 
sus redshift plane to characterize better the properties and evo- 
lution of typical AGNs and galaxies. The recent extension of 
the CDF-S survey from 2 Ms (Luo et al. 2008; hereafter LOS) 
to 4 Ms of exposure, via a large Director's Discretionary Time 
project, has now provided our most sensitive 0.5-8 keV view 
of the distant universe. These data, complemented by the re- 
cent w 3.3 Ms XMM-Newton observations in the CDF-S (Co- 
mastri et al. 2011), will enable detailed studies of AGN evo- 
lution, physics, and ecology as well as the X-ray properties of 
normal and starburst galaxies, groups and clusters of galaxies, 
large-scale structures, and Galactic stars. 

In this paper, we present Chandra source catalogs and data 
products derived from the full 4 Ms CDF-S data set as well as 
details of the observations, data reduction, and technical anal- 
ysis. We have made a number of methodological improve- 
ments in catalog production relative to past CDF catalogs. 
The structure of this paper is the following: in § 2 we de- 
scribe the observations and data reduction; in § 3 we detail the 
production of images, exposure maps, and the candidate-list 
catalog; in § 4 and § 5 we present the main and supplementary 
source catalogs as well as description of the adopted method- 
ology, respectively; in § 6 we perform simulations to assess 
the completeness and reliability of the main source catalog; 
in § 7 we estimate the background and sensitivity across the 
CDF-S and investigate the prospects for longer CDF-S expo- 
sures; and in § 8 we summarize the results of this work. 

Throughout this paper, we adopt a Galactic column density 
of A^H = 8.8 X 10''' cm-2 (e.g.. Stark et al. 1992) along the 
line of sight to the CDF-S. We use J2000.0 coordinates and 
a cosmology of //o = 70.4 km s"' Mpc"\ Om = 0.272, and 
f^A = 0.728 (e.g., Komatsu et al. 201 1). 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Observations and Observing Conditions 

Table 1 summarizes basic information for the 3 1 CDF-S ob- 
servations that were taken between 2010 March 18 and 2010 
July 22, which comprise the second 2 Ms exposure. The first 
2 Ms exposure consisted of 23 observations (see Table 1 of 
LOB for basic information) that were performed between 1999 



October 15 and 2007 November 4; the corresponding source 
catalogs were presented in LOS. 

All 54 CDF-S observations made use of the Advanced CCD 
Imaging Spectrometer imaging array (ACIS-I; Garmire et al. 
2003). The ACIS-I is comprised of four 1024 x 1024 pixel 
CCDs (CCDs 10-13; each has a pixel size of 0.492") and 
is optimized for imaging wide fields (with a field of view 
of 16.9' X 16.9' = 285.6 arcmin^). The focal-plane temper- 
ature was -110°C during the first two observations (1431-0 
and 1431-1; Giacconi et al. 2002; LOS) and -120°C during 
the others. The 10 early CDF-S observations between 1999 
November 23 and 2000 December 23 were taken in Faint 
mode (Giacconi et al. 2002; LOS); all the later CDF-S obser- 
vations as well as the earliest one (observation 1431-0) were 
taken in Very Faint mode in order to improve the screening of 
background events and thus increase the sensitivity of ACIS 
in detecting faint X-ray sources (Vikhlinin 2001). 

We inspected the background light curves for all 54 CDF-S 
observations using the Chandra Imaging and Plotting Sys- 
tem (ChlPS)^'* as well as EVENT BROWSER in the Tools for 
ACIS Real-time Analysis (TARA; Broos et al. 2000) software 
package. We find no significant flaring for all observations 
(the background is stable within w20% of typical quiescent 
Chandra values) except observation 1431-0, during which a 
mild flare with a factor of « 3 increase for « 5 ks occurred. 
We filtered the data on good-time intervals, removed the one 
mild flare, and obtained a total exposure time of 3.872 Ms for 
the 54 CDF-S observations. 

The entire CDF-S covers an area of 464.5 arcmin^; this is 
considerably larger than the AClS-1 field of view because the 
aim points and roll angles vary between observations. The 
average aim point, weighted using the 54 individual exposure 
times, is ajiooo.o = 03'>32'"28.06^ (Jjiooo.o = -27°48'26.4". 

2.2. Data Reduction 

Table 1 lists the versions of the Chandra X-ray Center 
(CXC) pipeline software used to process the basic archive 
data products for the 3 1 new observations (see Table 1 of LOS 
for the information for the first 23 observations). We closely 
followed LOS in reducing and analyzing the data and refer 
readers to LOS for details. Briefly, we utilized Chandra Inter- 
active Analysis of Observations (CIAO; we used CIAO 4.2 and 
CALDB 4.3.0) tools and custom software, including the tara 
package (version released on 2010 February 26), as appropri- 
ate. 

We reprocessed each level 1 observation with the CIAO tool 
ACiS_PROCESS_EVENTS to correct for the radiation damage 
sustained by the CCDs during the first few months of Chandra 
operations using a Charge Transfer Inefficiency (CTI) correc- 
tion procedure (Townsley et al. 2000, 2002)^'', to remove the 
standard pixel randomization which causes point spread func- 
tion (PSF) blurring, and to apply a modified bad-pixel file. 
We made use of a customized stripped-down bad-pixel file 
rather than the standard CXC bad-pixel file because the latter 
excludes « 6-7% of the ACIS-I pixels on which a large frac- 
tion of events are valid for source searching as well as pho- 
tometry and spectral analysis (see § 2.2 of LOS for details). 
Our bad-pixel screening removed w 1 .3% of all events. When 

See http://cxc.harvard.edU/ciao3.4/download/doc/chips_manual/ for the 
Chips reference manual. 

TARA is available at http://www.astro.psu.edu/xray/docs/TARA/. 
2* The CXC CTI correction procedure is only available for -120°C data 
and is thus not applied to observations 1431-0 and 1431-1. 
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TABLE 1 

Journal of New Chandra Deep Field-South Observations 



Obs. Start Exposure Aim Point'' Roll Angle*^ Pipeline 



Obs. ID 


(UT) 


Time" (ks) 


a (J2000.0) 


6 (J2000.0) 




V oHJll 




[The 23 observations made during the first 2 Ms exposure are listed in 


Table 1 of L08] 




12043 


2010 Mar 18, 01:39 


129.6 


03 32 28.78 


-27 48 52.1 


252.2 


8.2.1 


12123 


2010 Mar 21, 08:08 


24.8 


03 32 28.78 


-27 48 52.1 


252.2 


8.2.1 


12044 


2010 Mar 23, 11:31 


99.5 


03 32 28.55 


-27 48 5 1 .9 


246.2 


8.2.1 


12128 


2010 Mar 27, 13:08 


22.8 


03 32 28.55 


-27 48 5 1 .9 


246.2 


8.2.1 


12045 


2010 Mar 28, 16:38 


99.7 


03 32 28.32 


-27 48 5 1 .4 


240.2 


8.2.1 


12129 


2010 Apr 03, 15:21 


77.1 


03 32 28.33 


-27 48 51.4 


240.2 


8.2.1 


12135 


2010 Apr 06,09:36 


62.5 


03 32 28.01 


-27 48 50.2 


231.7 


8.2.1 


12046 


2010 Apr 08, 08:17 


78.0 


03 32 28.01 


-27 48 50.2 


231.7 


8.2.1 


12047 


2010 Apr 12, 13:21 


10.1 


03 32 27.80 


-27 48 48.9 


225.2 


8.2.1 


12137 


2010 Apr 16, 08:53 


92.8 


03 32 27.59 


-27 48 47.2 


219.2 


8.2.1 


12138 


2010 Apr 18, 12:40 


38.5 


03 32 27.59 


-27 48 47.3 


219.2 


8.2.1 


12055 


2010 May 15, 17:15 


80.7 


03 32 26.72 


-27 48 32.3 


181.4 


8.2.1 


12213 


2010 May 17, 14:22 


61.3 


03 32 26.69 


-27 48 31.1 


178.9 


8.2.1 


12048 


2010 May 23, 07:09 


138.1 


03 32 26.64 


-27 48 27.6 


171.9 


8.2.1 


12049 


2010 May 28, 18:58 


86.9 


03 32 26.61 


-27 48 24.4 


165.5 


8.2.1 


12050 


2010 Jun 03, 06:47 


29.7 


03 32 26.61 


-27 48 21.7 


160.2 


8.2.1 


12222 


2010 Jun 05, 02:47 


30.6 


03 32 26.61 


-27 48 21.7 


160.2 


8.2.1 


12219 


2010 Jun 06, 16:30 


33.7 


03 32 26.61 


-27 48 21.7 


160.2 


8.2.1 


12051 


2010 Jun 10, 11:30 


57.3 


03 32 26.63 


-27 48 19.2 


155.2 


8.2.1 


12218 


2010 Jun 11, 10:18 


88.0 


03 32 26.63 


-27 48 19.2 


155.2 


8.2.1 


12223 


2010 Jun 13, 00:57 


100.7 


03 32 26.63 


-27 48 19.2 


155.2 


8.2.1 


12052 


2010 Jun 15, 16:02 


110.4 


03 32 26.70 


-27 48 14.5 


145.7 


8.2.1 


12220 


2010 Jun 18, 12:55 


48.1 


03 32 26.70 


-27 48 14.5 


145.7 


8.2.1 


12053 


2010 Jul 05, 03:12 


68.1 


03 32 27.02 


-27 48 06.0 


127.0 


8.3 


12054 


2010 Jul 09, 11:35 


61.0 


03 32 27.02 


-27 48 06.1 


127.0 


8.3 


12230 


2010 Jul 11,03:52 


33.8 


03 32 27.02 


-27 48 06.0 


127.0 


8.3 


12231 


2010 Jul 12, 03:22 


24.7 


03 32 27.16 


-27 48 03.6 


121.2 


8.3 


12227 


2010 Jul 14, 21:04 


54.3 


03 32 27.16 


-27 48 03.7 


121.2 


8.3 


12233 


2010 Jul 16, 10:25 


35.6 


03 32 27.16 


-27 48 03.7 


121.2 


8.3 


12232 


2010 Jul 18, 19:53 


32.9 


03 32 27.16 


-27 48 03.7 


121.2 


8.3 


12234 


2010 Jul 22, 19:58 


49.1 


03 32 27.19 


-27 48 03.3 


120.2 


8.3 



Note. — The 4 Ms CDF-S consists of 54 observations, with the first 2 Ms exposure composed of 23 
observations (listed in Table 1 of LOS; not listed here to avoid repetition) and the second 2 Ms exposure 
composed of 31 observations (listed in this table; these 31 observations were all taken with the Very Faint 
mode). Right ascension has units of hours, minutes, and seconds, and declination has units of degrees, 
arcminutes, and arcseconds. 

^ Each of the 54 observations was continuous. We filtered the data on good-time intervals and removed one 
mild flare in observation 1431-0 (during the first 2 Ms exposure). The summed exposure time for the 54 

observations is 3.872 Ms. 

The aim points of the individual observations are the nominal ones taken from the Chandra archive. The 
average aim point, weighted by the 54 exposure times, is cijaooo.o = 03''32™28.06^ 5j2O0O.O - — 27°48'26.4". 

Roll angle, describing the orientation of the Chandra instruments on the sky, ranges from 0° to 360° and 
increases to the west of north (opposite to the sense of traditional position angle). 

The version of the CXC pipeline software used for the basic processing of the data. 



cleaning background events, we set CHECK_VF_PHA=YES in 
AClS_PROCESS_EVENTS for observations taken in Very Faint 
mode to utilize a 5 x 5 pixel event island to search for poten- 
tial cosmic-ray background events, which typically removes 
« 20-30% of the events of individual observations. 

We used the CIAO tool ACIS_DETECT_AFTERGL0W to re- 
move cosmic-ray afterglows, which is more stringent than 
the ciAO tool ACis_RUN_HOTPix that often fails to flag a 
substantial number of obvious cosmic-ray afterglows. Even 
ACiS_DETECT_AFTERGLOW fails to reject all afterglows. 
Working in CCD coordinates, we therefore utilized custom 
software to clean the data further by removing many addi- 
tional faint afterglows with 3 or more total counts occur- 
ring within 20 s (or equivalently 6 consecutive frames) on a 
pixel.'^^ We removed a total of 176 additional faint afterglows 

As shown later in Table 8, the fuU-band (i.e., 0.5-8 keV) mean back- 
ground rate of the 4 Ms CDF-S is 0.252 count Ms"' pixel"' , which translates 
into a coimt rate of 5.04 x 10"* counts per 20 s per pixel. Given such a low 
background count rate, the probabiUty of 3 or more counts (that are not asso- 
ciated with cosmic-ray afterglows) occurring within 20 s on a pixel by chance 



across the full 4 Ms dataset which, upon inspection, were 
isolated and not associated with apparent legitimate X-ray 
sources. 

As stated above, one significant deviation of our data re- 
duction from the CXC reduction of the CDF-S data set^^ is 
implementation of a customized stripped-down bad-pixel file, 
which retains an appreciable number of valid events (account- 
ing for w 5% of all events) that would have been discarded 
using the standard CXC bad-pixel file. As will be described 
in § 3.1 and § 3.2, when creating the CDF-S data products 
(e.g., merged X-ray images) and source catalogs (e.g., X-ray 
source positions), we registered individual X-ray observations 
to a common optical/radio astrometric frame (see § 3.1) and 
refined the absolute astrometry of the merged X-ray images 
and source positions using high-quality radio data (see § 3.2), 
thereby producing sharp merged X-ray images and accurate 

is negligible (2.54 x 10""). 

2^ The CXC CDF-S data products are available at 
http://cxc.harvard.edu/cda/Contrib/CDFS.html. 
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X-ray source positions (with < 0.2" astrometric shifts); in 
contrast, the CXC did not utilize muMwavelength data to reg- 
ister and refine X-ray astrometry. 

3. IMAGES, EXPOSURE MAPS, AND CANDIDATE-LIST CATALOG 

While following the general procedure described in § 3 
of LOS in the production of our source catalogs, we exten- 
sively made use of the ACIS Extract (AE; version released 
on 2010 February 26; Broos et al. 2010)^^ point-source anal- 
ysis software that appropriately computes source properties 
when multiple observations with different roll angles and/or 
aim points are being combined (such as those analyzed here). 
Significant improvements from the methodology of LOS in- 
clude, e.g., (1) utilization of AE polygonal source-count ex- 
traction regions that approximate the shape of the PSF and 
take into account the multi-observation nature of the data, 
and (2) utilization of a two-stage approach to source detec- 
tion, which filters candidate sources according to binomial 
no-source probabilities (i.e., probabilities of sources not being 
real considering their local backgrounds) calculated by AE. 

We first generated a candidate-list catalog of sources de- 
tected by WAVDETECT (Freeman et al. 2002) on the combined 
images (see § 3.2) at a false-positive probability threshold of 
10"^. We then pruned the candidate-list catalog to obtain a 
more conservative main catalog by removing low-significance 
source candidates, according to the AE-computed binomial 
no-source probabilities. As detailed later in § 3.2 and § 4, 
this approach not only produces source catalogs that are of 
similar quality to those produced by running WAVDETECT at 
the more typical false-positive probability threshold of 10~* 
or 10"^ used in previous CDF studies (e.g., Alexander et al. 
2003, hereafter A03; Lehmer et al. 2005, hereafter LOS; L08), 
but also allows for flexibility in including additional legiti- 
mate sources that fall below the lO"*" or 10"^ threshold. This 
procedure has previously been employed in similar forms in 
a number of studies (e.g., Getman et al. 2005; Nandra et al. 
2005; Laird et al. 2009; Lehmer et al. 2009). 

3.L Image and Exposure Map Creation 

To construct the combined event file we initially ran 
WAVDETECT at a false-positive probability threshold of 10"^ 
on the individual cleaned 0.5-8 keV image of each obser- 
vation to generate initial source lists and used AE to deter- 
mine centroid positions of each detected source. We then 
registered the observations to a common astrometric frame 
by matching X-ray centroid positions to optical sources de- 
tected in deep /?-band images taken with the Wide Field 
Imager (WFT) mounted on the 2.2-m Max Planck Gesell- 
shaft/European Southern Observatory (ESO) telescope at La 
Silla (see § 2 of Giavalisco et al. 2004). We have manu- 
ally shifted all the WFI /?-band source positions by 0.175" 
in right ascension and -0.284" in declination (also see Luo 
et al. 2010; hereafter LIO) to remove the systematic off- 
sets between the optical positions and the radio positions of 
sources in the Very Large Array (VLA) 1.4 GHz radio cat- 
alog presented in Miller et al. (2008).^*^ We did not di- 
rectly match X-ray centroid positions to the VLA radio cat- 
alog because, for some observations, there are too few com- 

^' See http://www.astro.psu.edu/xray/docs/TARA/ae_users_guide.html 
for details on ACIS Extract. 

Throughout this paper, we used the Scr VLA 1 .4 GHz radio catalog (N. 
A. Miller 2010, private communication) that has a hmiting flux density of 
« 40 fiiy. 



mon sources between the X-ray and radio source lists to en- 
sure a robust astrometric solution, owing to the relatively low 
radio-source density and relatively small numbers of Chan- 
dra sources detected in individual observations. However, as 
detailed in § 3.2, we are able to lock the absolute astrometry 
of the combined X-ray images to the VLA radio catalog be- 
cause of the larger number of X-ray sources detected. We per- 
formed X-ray/^-band matching and astrometric re-projection 
using the CIAO tools reproject_aspect and wcs_update 
with a 3" matching radius and a residual rejection limit^' of 
0.6". Typically, 60-150 X-ray//?-band matches were used 
in each observation for the astrometric solutions. When us- 
ing WCS_UPDATE, linear translations range from 0.032" to 
0.525", rotations range from -0.048° to 0.035°, and scale 
changes range from 1.00004 to 1.00145. Individual regis- 
trations are accurate to « 0.3". We then reprojected all the 
observations to the frame of observation 2406, which is one 
of the observations that requires the smallest translation to be 
aligned with the optical astrometric frame; however, we note 
that it does not matter which observation is used as the refer- 
ence frame for reprojection once each observation is analyzed 
consistently. 

We utilized the CiAO tool dmmerge to produce a merged 
event file by combining the individual event files. We con- 
structed images from this merged event file using the standard 
ASCA grade set (ASCA grades 0, 2, 3, 4, 6) for three standard 
bands: 0.5-8.0 keV (full band; FB), 0.5-2.0 keV (soft band; 
SB), and 2-8 keV (hard band; HB).^^ pigure 1 shows the raw 
full-band image. We generated effective-exposure maps for 
the three standard bands following the basic procedure out- 
lined in §3.2 of Hornschemeier et al. (2001) and normalized 
them to the effective exposures of a pixel located at the av- 
erage aim point. This procedure takes into account the ef- 
fects of vignetting, gaps between the CCDs, bad-column fil- 
tering, bad-pixel filtering, and the spatial and time dependent 
degradation in quantum efficiency due to contamination on 
the ACIS optical-blocking filters; thus, the derived effective 
exposures are typically smaller than the nominal exposures 
(i.e., durations of observations). When creating the effective- 
exposure maps, we assumed a photon index of F = 1 .4, the 
slope of the cosmic 2-10 keV X-ray background (e.g., Mar- 
shaUetal. 1980; Gendreau et al. 1995; Hasinger et al. 1998; 
Hickox & Markevitch 2006). Figure 2 shows the full-band 
effective-exposure map, and Figure 3 displays the survey solid 
angle as a function of the minimum full-band effective expo- 
sure. According to Fig. 3, about 52% and 38% of the CDF-S 
field has a full-band effective exposure greater than 2 Ms 
and 3 Ms, respectively; the maximum effective exposure is 
3.811 Ms, which is slightly smaller than the 3.872 Ms total 
exposure since the locations of the aim points of individual 
observations vary. For a given full-band effective exposure, 
the survey solid angle is up to a factor of w 1.5 times larger 
than that of the 2 Ms CDF-S (L08; Fig. 3, dash-dot curve) at 
the low end of effective exposure (< 1.5 Ms), and it is much 

^' This is a parameter used in WCS_UPDATEto remove source pairs based 
on pair positional offsets. 

32 We compared X-ray source catalogs made with the two upper energy 
cuts of 7 keV and 8 keV (i.e., the set of energy bands of 0.5-7.0, 0.5-2.0, 
and 2-7 keV versus the set of energy bands adopted here; see § 3.2 and § 4.1 
for the details of catalog production). We foimd no clear statistical difference 
between catalogs; the X-ray sources that are unique in each catalog are faint 
(i.e., close to or right on source-detection limits) and account for only ^ 3% 
of all detected sources. We thus adopted the traditional standard bands (i.e., 
using the upper energy cut of 8 keV) to maintain continuity with past catalogs 
(e.g., A03; LOS; LOS). 
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larger than that of the 2 Ms CDF-S above 1.5 Ms effective 
exposure. Thus in addition to the fact that the 4 Ms CDF-S 
can detect new sources that have lower fluxes than the 2 Ms 
sources, it can also detect new sources that have a similar flux 
distribution to the 2 Ms sources over as much as 50% more 
area. 

We followed § 3.3 of Baganoff et al. (2003) to con- 
struct exposure-corrected smoothed images. We first pro- 
duced the raw images and effective-exposure maps in the 
0.5-2.0 keV, 2-A keV, and 4-8 keV bands, using the afore- 
mentioned procedures. We then adaptively smoothed the 
raw images and effective-exposure maps using the CIAO 
tool CSMOOTH (EbeHng, White, & Rangarajan 2006). Fi- 
nally, we divided the smoothed images by their correspond- 
ing smoothed effective-exposure maps and combined the 
exposure-corrected smoothed images together to produce a 
full-band color composite, as shown in Fig. 4 (note that this 
color composite is not background- subtracted); an expanded 
view of the central 8' x 8' region is also shown in Fig. 4. Note 
that we ran wavdetect only on the raw images for source 
searching, although many detected X-ray sources appear more 
clearly in the adaptively smoothed images. 

3.2. Candidate-List Catalog Production 

We ran wavdetect on each combined raw image in the 
three standard bands^^ to perform source searching and to 
construct a candidate-list catalog, using a "\/2 sequence" of 
wavelet scales (i.e., 1, Vl, 2, 2^2, 4, 4V^, 8, 8^2, and 
16 pixels) and a false-positive probability threshold of 10"^. 
We expect the use of a false-positive probability threshold of 
10"^ to introduce a non-negligible number of spurious sources 
that have <2-3 source counts. However, as pointed out 
by Alexander et al. (2001), using a more stringent source- 
detection threshold (e.g., 10~^, 10"^, or 10"^) can lose an ap- 
preciable number of real sources. In § 4, we create a more 
conservative main catalog by determining the detection sig- 
nificances of each candidate-list source in the three standard 
bands and discarding sources with significances below an 
adopted threshold value. 

Our candidate-list catalog consists of 892 X-ray source can- 
didates; each candidate was detected in at least one of the 
three standard bands with WAVDETECT at a false-positive 
probabiUty threshold of 10~^. We adopted, in order of pri- 
ority, full-band, soft-band, or hard-band source positions for 
candidate sources. We performed cross-band matching using 
a 2.5" matching radius for sources within 6' of the average 
aim point (i.e., off-axis angle 6 < 6') and a 4.0" matching ra- 
dius for sources located at larger off-axis angles (i.e., 9 > 6'). 
The choice of these matching radii was made based on in- 
spection of histograms that show the number of matches as 
a function of angular separation (e.g., see §2 of Boiler et al. 
1998). With these matching radii, the mismatch probability 
is ~ 1% over the entire field. We removed a few duplicate 
sources due to false matches near the edge of the field through 
visual inspection.^"* 

We note that WAVDETECT was run on the combined raw images where 
the average aim point (given in § 2) is a good approximation of the image cen- 
ter for the purpose of computing PSFs. Given that we used multiple wavelet 
scales, the Mexican-Hat wavelet patterns (adopted by wavdetect) provide 
reasonable first-order approximations of the multi-observation PSFs. 

^ For a few sources that lie near the edge of the field, the offset between 
the X-ray positions determined from different bands by wavdetect is > 4" ; 
such a source will be counted twice (i.e., treated as two sources) according to 



We improved the above WAVDETECT source positions uti- 
lizing the centroid and matched-filter positions computed by 
AE. The matched-filter positions are obtained by correlating 
the full-band image in the neighborhood around each source 
with the source's combined PSF. The combined PSF is gener- 
ated by combining the individual PSFs of a source for each 
relevant observation, weighted by the number of detected 
counts. This technique takes into account the fact that, due to 
the complex PSF at large off-axis angles, the X-ray source po- 
sition is not always located at the peak of the X-ray emission. 
The WAVDETECT, centroid, and matched-filter positions have 
comparable accuracy on-axis, while the matched-filter posi- 
tions have better accuracy off-axis. Thus, we adopted centroid 
positions for sources with 9 <S' and matched-filter positions 
for sources with 9 > 8'. 

We refined the absolute astrometry of the raw X-ray im- 
ages by matching the candidate-list sources to the 5cr VLA 
1.4 GHz radio-catalog sources (see § 3.1). There are 359 ra- 
dio sources across the CDF-S field with positions accurate 
to < 0.1". We performed cross-matching between the 892 
candidate-list catalog X-ray sources and the 359 radio sources 
in the field using a 2" matching radius and found 141 matches. 
We estimated the expected false matches by manually shifting 
the X-ray source positions in right ascension and declination 
by ±(5-60") in steps of 5" (i.e., in unique directions) and 
recorrelating with the radio sources. The average number of 
false matches is « 2.3 (« 1 .7%) and the median offset of these 
false matches is 1.41". Of the 141 matches, we identified five 
extended radio sources upon inspecting the radio image. We 
excluded two of these five extended radio sources for the as- 
trometry refinement analysis because these two matches are 
spurious with positional offsets greater than 1 .5" (see § 4.2 
for more details on these two extended radio sources); the 
other three matches are robust with small positional offsets 
(< 0.7") and were included for the subsequent analysis. Us- 
ing these 139 matches, we found small shift and plate-scale 
corrections when comparing the X-ray and radio source posi- 
tions and applied these corrections to all the combined X-ray 
images and source positions, which results in small (< 0.2") 
astrometric shifts. 

We utilized AE to perform photometry for the candidate-list 
catalog sources. Compared to "traditional" circular-aperture 
photometry (e.g., L08), the most important difference in the 
AE-computed photometry is the use of polygonal source- 
extraction regions.^^ AE models the Chandra High Reso- 
lution Mirror Assembly (HRMA) using the MARX^^ ray- 
tracing simulator (version 4.4.0) to obtain the PSF model. 
It then constructs a polygonal extraction region that ap- 
proximates the « 90% encircled-energy fraction (FEE) con- 
tour of a local PSF measured at 1.497 keV (note that AE 
also constructs PSFs at energies of 0.277, 4.510, 6.400, and 
8.600 keV). When dealing with crowded sources having over- 
lapping polygonal extraction regions, AE utilized smaller ex- 
traction regions (corresponding to R:!40-75% EEFs) that were 
chosen to be as large as possible without overlapping. Less 

our matching approach (i.e., a 4.0" matching radius at > 6'). We removed 
the duplicated sources in these few cases. 

The polygonal source-extraction regions typically become more non- 
circular toward larger off-axis angles. In particular, the source-extraction re- 
gions for crowded sources at large off-axis angles are reduced from 90% 
to ^ 40-75% encircled-energy fractions (EEFs) and thus represent the most 
dramatic examples of deviation from circular apertiffes (see, e.g.. Fig. 6 of 
Broos et al. 2010 for such an example). 

3* MARX is available at http://space.mit.edu/CXC/MARX/index.html. 
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Fig. 1 . — Full-band (0.5-8.0 keV) raw image of the 4 Ms CDF-S displayed with Unear gray scales. The segmented boimdary surrounding the image shows the 
coverage of the entire CDF-S. The large polygon, the rectangle, and the central small polygon indicate the regions for the GOODS-S (Giavalisco et al. 2004), the 
planned CANDELS GOODS-S (5-orbit HST/WVCS; see § 8 for more details about CANDELS), and the Hubble Ultra Deep Field (UDF; Beckwith et al. 2006), 
respectively. The central plus sign indicates the average aim point, weighted by exposure time (see Table 1). The pale ring-like area near the field center is caused 
by the ACIS-I CCD gaps in which the effective exposures are lower than in the nearby non-gap areas (see Fig. 2). The apparent scarcity of soitfces near the field 
center is mainly due to the small size of the on-axis PSF (see Figs. 4 and 13 for clarification). 



than 6% of the 892 candidate-Hst sources are crowded by this 
definition. For background extraction, we adopted the AE 
"BETTER_BACKGROUNDS" algorithm. This algorithm 
models the spatial distributions of flux for the source of inter- 
est and its neighboring sources using unmasked data. It then 
computes local background counts within background regions 
that subtract contributions from the source and its neighbor- 
ing sources. In our AE usage, the background-extraction 
region is typically a factor of w 16 larger than the source- 
extraction region and contains at least 1 00 background counts. 
As discussed in § 7.15 of the AE manual, AE also imposes 
an expUcit requirement that the uncertainty in the estimate 
of net counts be dominated by the uncertainty in the ex- 



tracted source counts in order to ensure photometric accu- 
racy; this requirement leads to enlargement of background re- 
gions/counts when necessary. As a result, the median number 
of full-band background counts extracted for the main-catalog 
sources (see § 4.4) is 780, with an interquartile range of 
278-2621. This algorithm produces accurate background ex- 
tractions, which are particularly critical for crowded sources. 
For sources that are not crowded, this algorithm produces es- 
sentially the same background-extraction results as the tra- 
ditional AE "EXTRACT_BACKGROUNDS" algorithm; the 
latter algorithm computes local background counts by mask- 
ing all the sources and then searching around each source for 
the smallest circular region that contains a desired number of 
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Fig. 2. — Full-band (0.5-8.0 keV) effective-exposure map of the 4 Ms 
CDF-S displayed with linear gray scales that are indicated by the inset scale 
bar (effective exposure times are in units of seconds). The darkest areas rep- 
resent the highest effective exposure times, with a maximum of 3.811 Ms. 
The distributions of the ACIS-I CCD gaps can be clearly identified (indicated 
by the radial trails). The regions and the plus sign are the same as those in 
Fig. 1. 
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Fig. 3. — Plot of survey solid angle as a function of minimum fuU-band 
effective exposure for the 4 Ms CDF-S (solid curve). The maximum exposure 
is 3.811 Ms. The vertical dotted line indicates an effective exposure of 2 Ms. 
Approximately 242.3 arcmin^ (Ri 52%) of the CDF-S survey area has > 2 Ms 
effective exposure. For comparison, the 1 Ms CDF-S result {dashed curve) 
and the 2 Ms CDF-S result (dash-dot curve), both of which are obtained using 
the procedures in this paper, are also shown in the plot. 

background counts. AE analyzes individual observations in- 
dependently (including, e.g., the use of MARX for PSF mod- 
eling and source and background extractions) and merges the 
data to produce photometry for each source. The resulting 
combined PSFs at 1.497 keV have typical FWHMs of 0.68", 
1.07", 1.76", 2.79", and 3.61" at off-axis angles of 1', 3', 5', 
7', and 9', respectively; these FWHM values represent typical 

For this work, we did not use the optional AE 
"MERGE_FOR_PHOTOMETRY" algorithm, as discussed in Broos 
et al. (2010), that allows AE to discard some extractions during a merge of 
AE products from individual observations. 



angular resolutions of the 4 Ms images. 

AE estimates an energy-dependent aperture correction for 
each source and apphes the correction to the effective area 
calibration file used for spectral modeling. For this work, 
we chose to apply aperture corrections to the background- 
subtracted photometry as follows. For the soft (hard) band, we 
derived an effective PSF fraction for each source by weight- 
ing PSF measurements at 1.497 (4.510) keV by the exposures 
for the individual observations. Given that the full band is 
a combination of the soft and hard bands, we derived the 
full-band effective PSF fraction based on the derived soft- 
band and hard-band effective PSF fractions: (1) if a source 
was detected both in the soft and hard bands, we derived 
the full-band effective PSF fraction by weighting the soft and 
hard-band effective PSF fractions with the soft and hard-band 
background- subtracted counts; (2) if a source was detected 
in the soft or hard band (but not both), we set the full-band 
effective PSF fraction to the soft or hard band effective PSF 
fraction, respectively; and (3) if a source was detected in nei- 
ther the soft band nor the hard band, we took the average of 
the soft and hard band effective PSF fractions as the full-band 
effective PSF fraction. The median aperture corrections for 
the full, soft, and hard bands are 0.875, 0.898, and 0.826, re- 
spectively. We then applied aperture corrections by dividing 
the background-subtracted source counts by the derived effec- 
tive PSF fractions. Since our candidate-list catalog was con- 
structed using WAVDETECT with a liberal false-positive prob- 
ability threshold of 10"^, many candidate sources have < 2-3 
(background-subtracted) source counts. In the next section, 
we evaluate the reliability of candidate sources on a source- 
by-source basis to produce a more robust main source catalog . 

4. MAIN CHANDRA SOURCE CATALOG 
4.1. Selection of Main-Catalog Sources 

As discussed above, we expect our candidate-list catalog 
of 892 X-ray sources to include a significant number of false 
sources since we ran WAVDETECT at a liberal false-positive 
probability threshold of 10"^. If we conservatively treat the 
three standard-band images as independent, we can estimate 
the number of expected false sources in the candidate-hst cat- 
alog for the case of a uniform background by multiplying the 
WAVDETECT threshold of 10"^ by the sum of pixels in the 
three bands (i.e., w 2.07 x 10^). However, such a false-source 
estimate is conservative, since over the majority of the field, 
a single pixel will not be considered a source-detection cell. 
In particular, at large off-axis angles wavdetect suppresses 
fluctuations on scales smaller than the PSF. As quantified in 
§ 3.4.1 of A03, the number of false-sources is likely fs2-3 
times smaller than the above conservative estimate. We refer 
readers to § 6.2 for relevant discussions. 

To produce a more reliable main Chandra source catalog, 
we evaluated for each source the binomial probability P that 
no source exists given the measurements of the source and lo- 
cal background. As discussed in § 5.10.3 of the AE manual 
(also see Appendix A2 of Weisskopf et al. 2007 for further de- 
tails), the binomial no-source probabiUty P can be calculated 
using the following equation: 

N , 

Pix > 5) = V '■ — p^(i -pf~^. (1) 

In this equation, S is the total number of counts in the 
source-extraction region without subtraction of the back- 
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Fig. 4. — (a) Chandra "false-color" image of the 4 Ms CDF-S, which is a color composite of the exposure-corrected and adaptively smoothed images in the 
0.5-2.0 keV (red), 2—4 keV (green), and 4—8 keV (blue) bands, (b) An expanded view of the Chandra "false-color" image of the central 8' X 8' region (note that 
a slightly different contrast ratio from that for the full image is used here in order to render the faint sources more clearly). The apparent smaller size and lower 
brightness of sources near the field center is due to the smaller size of the on-axis PSF. The regions and the plus sign are the same as those in Fig. 1 . 
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ground counts Bsic in this region; A' = 5 + where Bext is 
the total extracted background counts within a background- 
extraction region that is typically a factor of « 16 larger than 
the source-extraction region in our AE usage (see § 3.2); 
and /> = l/{l+BACKSCAL) is the probability that a photon 
lies in the source-extraction region (thus contributing to 5), 
where BACKSCAL = -Sext/^src with a typical value of s::^ 16, 
as stated earUer. P is computed by AE in each of the three 
standard bands. For a source to be included in our main cat- 
alog, we required P < 0.004 in at least one of the three stan- 
dard bands. We identified multiwavelength counterparts for 
the X-ray sources (see § 4.3) and studied the identification 
rate as a function of the P value, given that X-ray sources 
without identifications in ultradeep multiwavelength data are 
more Ukely to be false detections (see, e.g., LIO). The re- 
quirement of P < 0.004 was empirically chosen as a com- 
promise to keep the fraction of potential false sources small 
while recovering the largest number of real sources. Using 
this criterion of P < 0.004, our main catalog contains a total 
of 740 sources. We note that for a different choice of source- 
detection criterion of P < 0.01, a total of 33 additional sources 
with 0.004 < P < 0.01 would be included; however, only 
« 64% (i.e., 21) of these 33 sources have multiwavelength 
counterparts, as opposed to an identification rate of « 97% 
for the main catalog (see § 4.3). We refer readers to § 6.2 
for a detailed discussion on the completeness and reliability 
of the main catalog based on simulations. 

Our adopted cataloging procedure, with the utilization 
of AE, has a number of advantages over a "traditional" 
WAVDETECT-only approach: (1) the more detailed treatment 
of complex source-extraction regions (i.e., using polygonal 
regions, as opposed to elhptical apertures, to simulate the 
PSF) that is more suitable for the case of multiple observa- 
tions with different aim points and roll angles, (2) the better 
source-position determination that maximizes the signal-to- 
noise ratio and leads to more accurate count estimates, (3) the 
more careful background estimates that take into account the 
effects of all the neighboring sources and CCD gaps, and (4) 
the more immediately transparent mathematical criterion (i.e., 
the binomial probability) that is utilized for source detection. 
We will demonstrate below that our adopted procedure recov- 
ers almost all of the sources detected with WAVDETECT at a 
false-positive probabiUty threshold of 10"^ and a significant 
number of additional real sources detected at 10"^. 

In order to give a more detailed WAVDETECT-based per- 
spective on source significance, we also ran wavdetect 
on the three standard-band images at false-positive proba- 
bility thresholds of lO"*", 10"^, and 10"**, and found detec- 
tions for 659 (73.9%), 569 (63.8%), and 502 (56.3%) of 
the 892 candidate-Ust catalog sources, respectively. Among 
the 152 candidate-list sources that failed the selection cut of 
P < 0.004, and thus were not included in the main catalog, 40 
(« 4.5% of the 892 candidate-list sources) had WAVDETECT 
false-positive probability detection thresholds of < 10"^. 
Meanwhile, our main catalog includes 121 sources that had 
minimum WAVDETECT probabiUties of 10~^.^^ Therefore, our 
adopted procedure, as opposed to a direct WAVDETECT-based 

The minimum WAVDETECT probability represents the WAVDETECT sig- 
nificance of a source, witli lower values indicating higher significances. For 
example, if a source was detected with WAVDETECT in at least one of the 
three standard bands at a false-positive probability threshold of Itr^ but was 
not detected in any of the three standard bands at a threshold of 10"**, then 
the minimum wavdetect probability of this source is 10~^. 



approach, has a "net gain" of 8 1 sources. We note that a larger 
net gain of sources could be achieved if we adopted a less 
conservative no-source probability cut (e.g., P < 0.01) at the 
expense of introducing more spurious sources. 

Figure 5 shows the fraction of candidate-list sources in- 
cluded in the main catalog and the 1 - P distribution of 
candidate-list sources as a function of the minimum WAVDE- 
TECT probability. The fraction of candidate-list sources in- 
cluded in the main catalog is 98.2%, 85.1%, 76.7%, and 
51.9% for a minimum WAVDETECT probabihty of 10"^, 10"^, 
10"^, and 10"^, respectively. As shown later in § 4.3, we 
find that 716 (96.8%) of the 740 main-catalog sources have 
secure multiwavelength counterparts (with a false-matching 
probability of « 2. 1 %), where the identification rate is 98.1% 
(90.1%) for the 619 (121) sources with a minimum WAVDE- 
TECT probability of < 10"^ (10"^) in the main catalog. Given 
the relatively small false-matching rate, the above high iden- 
tification rates indicate that the vast majority of the main- 
catalog sources are real X-ray sources (see, e.g., LIO). Thus, 
our main-catalog selection provides an effective identification 
of real X-ray sources including those falling below the tradi- 
tional 10"^ WAVDETECT searching threshold. 

4.2. X-ray Source Positional Uncertainty 

As in § 3.2, we cross matched the 740 main-catalog sources 
with the 359 radio sources in the field using a 2" matching ra- 
dius and found 135 matches. We estimated on average 2.0 
(w 1 .5%) false matches and a median offset of 1 .45" for these 
false matches. Figure 6(a) shows the positional offset be- 
tween the X-ray sources and their radio counterparts as a func- 
tion of off-axis angle. The median positional offset is 0.24". 
There are three sources in Fig. 6(a) that have positional offsets 
greater than 1.5": (1) the one with the largest offset (1.97") 
mistakenly matches to one of the two lobes of a radio galaxy 
due to the fact that the radio core, which is likely the real 
counterpart for this X-ray source, was not detected in the radio 
catalog; (2) the one with the second largest offset (1.88") is 
likely a false match because such an offset is much larger than 
its expected positional uncertainty [see eq. (2) below] con- 
sidering its off-axis angle (11 .6') and source-counts (« 100); 
and (3) the one with the third largest offset (1.62") has one 
radio source, which is the core of a radio galaxy, and a few 
optical sources within its 2" radius, with the radio counter- 
part not matching to the Ukely real optical counterpart of this 
X-ray source (thus being a false match). Excluding the above 
three sources, we then estimated X-ray positional uncertain- 
ties using the remaining 132 X-ray detected radio sources. 
Figure 6(b) shows the positional residuals between the X-ray 
and radio positions for these 132 sources; the "scatter cloud" 
of positional residual appears circular, with no residual dis- 
tortions. As shown in Fig. 6(a), there are clear off-axis angle 
and source-count dependencies for these 132 sources, with 
the former due to the degradation of the Chandra PSF at large 
off-axis angles and the latter due to statistical limitations in 
finding the centroid of a faint X-ray source. Implementing 
the parametrization provided by Kim et al. (2007),""^ we de- 
rived an empirical relation for the positional uncertainty of 

3' We note that 6 (i.e., 141 - 135 = 6; also see § 3.2) candidate-Hst X-ray 
sources that have a radio counterpart were not included in the main catalog; 
these sources are likely real X-ray sources that fail to satisfy our relatively 
stringent source-selection criterion of P < 0.004 (see § 4.1). 

We note that the Kim et al. (2007) parametrization fits our data ade- 
quately (i.e., the AE-derived positions and photometry), although it was orig- 
inally based on WAVDETECT-derived positions and photometry. 
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Fig. 5. — The fraction of sources in the candidate-list catalog with an AE binomial no-source probability P < 0.004, which were included in the main 
catalog, as a function of minimum WAVDETECT probability'^ (shown as five-pointed stars). The number of sources with P < 0.004 versus the number of 
candidate-list catalog sources detected at each minimum WAVDETECT probability are annotated in the figure (note that, in this figure, 502+67+90+233=892 and 
493+57+69+121=740). The fraction of candidate-list catalog sources included in the main catalog falls from 98.2% to 51.9% between minimum WAVDETECT 
probabilities of 10"* and 10"'. Shown in the insets are the histograms of 1 - P for the candidate-list catalog sources at each minimum WAVDETECT probability, 
with shaded areas highlighting those included in the main catalog (i.e., having l-P > 0.996). 



our X-ray sources by fitting to these 132 X-ray sources that 
have radio counterparts within a radius of 1 .5". The relation 
is 

log Ax = 0.04846i-0.43561ogC-i-0.1258, (2) 

where Ax is the X-ray positional uncertainty in arcseconds, 9 
is the off-axis angle in arcminutes, and C is the source counts 
in the energy band where the source position was determined 
(see the description of Columns 8-16 of the main catalog in 
§ 4.4 for details on photometry calculation). We set an up- 
per limit of 2000 on C since the positional accuracy does not 
improve significantly above that level. As a guide to the de- 
rived relation, we show positional uncertainties for C = 20, 
200, and 2000 in Fig. 6. The stated positional uncertainties 
are for the w 68% confidence level, which are smaller than 
the WAVDETECT positional uncertainties, particularly at large 
off-axis angles, due to our adopted positional refinement. In 
Figure 7, we show the distributions of positional offset in four 
bins of X-ray positional uncertainty, as well as the expected 
false matches assuming a uniform spatial distribution of radio 
sources. For each histogram in Fig. 7, as expected, > 65% of 
the positional offsets between the X-ray sources and their ra- 
dio counterparts are less than the corresponding median X-ray 
positional uncertainty. 

Owing to the factor of « 2 increase in exposure/source 
counts from 2 Ms to 4 Ms, the areas of source positional error 



regions are expected to be reduced by w 30% on average (see 
§ 4.2 of LIO). We thus compared our positional uncertain- 
ties with the positional uncertainties for the 440 main-catalog 
sources that were previously detected in the L08 main cata- 
log (see the description of Column 59 of the main catalog in 
§ 4.4). We find a median ratio of 0.82 between our and the 
L08 positional uncertainties^' (corresponding to a median ra- 
tio of 0.67 between areas of our and the L08 positional error 
regions); such an improvement is in agreement with the above 
expectation. We also cross matched the 462 L08 main-catalog 
sources with the 359 radio sources in the field using a match- 
ing radius of 2", taking into account the systematic positional 
offsets between the optical catalogs and the VLA radio cata- 
log (see § 3.1). The median positional offset is 0.40" between 
the L08 main-catalog sources and their radio counterparts for 
a total of 94 matches, as opposed to 0.24" in our case. This 
significant improvement is not only because of the improved 
photon statistics, but also because we locked the astrometry of 
the combined X-ray images to the VLA radio sources rather 
than the WFI /?-band sources that were adopted by L08. 



'^^ In L08, the X-ray positional uncertainties are quoted at the Ri 85% con- 
fidence level. For straightforward comparison, we thus adopted the ^ 68% 
confidence-level positional uncertainties reported in Table 2 of LIO that were 
used in the LIO likelihood-ratio matching procedure (see § 4.3 for more de- 
tails). 
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Fig. 6. — (a) Positional offset vs. off-axis angle for the 135 main-catalog 
sources that have counterparts in the 5<j VLA 1.4 GHz radio catalog using 
a matching radius of 2" (see § 4.2 for descriptions of the three sources with 
> 1.5" positional offsets). Red filled, green filled, blue filled, and black open 
circles represent X-ray sources with > 2000, > 200, > 20, and < 20 counts 
in the energy band where the source position was determined, respectively. 
The red dotted curve shows the running median of positional offset in bins 
of 2'. The horizontal dashed line indicates the median offset (1.45") of the 
expected false matches. We used these data to derive the 68% confidence- 
level X-ray source positional uncertainties, i.e., eq. (2). Three solid curves 
indicate the 68% confidence-level positional uncertainties for sources with 
20, 200 and 2000 counts, (b) Positional residuals between the X-ray and 
radio positions for the 132 main-catalog sources that have radio counterparts 
within a radius of 1.5" [see Panel (a)]. Red and black filled circles indicate 
sources with an off-axis angle of < 6' and > 6', respectively. A large blue 
circle with a radius of 0.5" is drawn at the center as a guide to the eyes, [see 
the electronic edition of the Supplement far a color version of this figure.^ 
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Fig. 7. — Histograms showing the distributions of positional offset for 
the 135 main-catalog sources that have counterparts in the 5ct VLA 1.4 GHz 
radio catalog using a matching radius of 2". These 135 sources were di- 
vided into four bins according to their positional uncertainties estimated using 
eq. (2): 0"-0.25", 0.25"-0.50", 0.50"-0.75", and 0.75"-l". The vertical 
dashed line in each panel indicates the median X-ray positional uncertainty 
in each bin. The dotted line shows the total expected number of random radio 
sources as a function of the positional offset. < 35% of the radio counterpaits 
lie beyond the median X-ray positional uncertainty in each bin. 

4.3. Multiwavelength Identifications 

We utilized the likelihood-ratio matching procedure pre- 
sented in § 2 of LIO to identify the optical/near- 
infrared/infrared/radio (ONIR) counterparts for the main- 
catalog X-ray sources. Briefly, the likelihood-ratio technique 
(e.g., Sutherland & Saunders 1992; Ciliegi et al. 2003; Brusa 
et al. 2005, 2007) searches for probable counterparts taking 
into account the positional accuracy of both the ONIR and 
Chandra X-ray sources and also the expected magnitude dis- 
tribution of the counterparts. Compared to a simple match- 
ing method that searches for the nearest counterpart within 
a given radius, the likelihood-ratio method significantly re- 
duces the false-match probability toward faint ONIR magni- 
tudes (see, e.g., § 2.4 of LIO). 

We used seven ONIR catalogs for identification purposes 
(see Table 1 of LIO for further details): 

1. The ESO 2.2-m WFl /?-band catalog (denoted as 
"WFl"; Giavalisco et al. 2004), with a Sa Hmiting AB 
magnitude (Oke & Gunn 1983) of 27.3; 

2. The GOODS-S Hubble Space Telescope (HST) ver- 
sion r2.0z z-band catalog (denoted as "GOODS-S"; Gi- 
avalisco et al. 2004), with a 5a limiting AB magnitude 
of 28.2; 

3. The GEMS HST z-hmd catalog (denoted as "GEMS"; 
Caldwell et al. 2008), with a 5a limiting AB magnitude 
of 27.3; 

4. The GOODS-S MUSIC catalog (denoted as "MUSIC"; 
Grazian et al. 2006; we used the ^T-selected sources 
in the V2 catalog that was presented in Santini et al. 
2009) based on the Retzlaff et al. (2010) VLT/ISAAC 
data, with a limiting /T-band AB magnitude of 23.8 (at 
90% completeness); 

5. The MUSYC /T-band catalog (denoted as "MUSYC"; 
Taylor et al. 2009), with a 5cr limiting AB magnitude 
of 22.4; 
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6. The SIMPLE Spitzer/IRAC 3.6 /xm catalog (denoted as 
"SIMPLE"; Damen et al. 2011), with a 5a Umiting AB 
magnitude of 23.8; and 

7. The VLA 1.4 GHz radio catalog (denoted as "VLA"; 
Miller et al. 2008), with a 5a limiting flux density of 
w 40 ^Jy. 

As mentioned in § 3.1, we find systematic positional off- 
sets between the optical/near-infrared catalogs and the radio 
catalog and have chosen to shift all the optical and near- 
infrared/infrared source positions throughout this paper by 
0.175" in right ascension and -0.284" in declination to be 
consistent with the radio astrometry. 

We found that 716 (96.8%) of the 740 main-catalog sources 
have ONIR counterparts. For an X-ray source having multi- 
ple counterparts from the likelihood-ratio matching (108 such 
cases), we chose a primary counterpart from, in order of pri- 
ority, file VLA, GOODS-S, GEMS, MUSIC, WFI, MUSYC, 
or SIMPLE catalog. This order is chosen based on several 
related factors: the positional accuracy, angular resolution 
(to minimize any blending effects), false-match probabiHty, 
and catalog depth. Manual adjustments were made to a few 
sources based on visual inspection (e.g., we selected the op- 
tical position rather than the VLA radio position if the radio 
counterpart is clearly extended; see § 2.3 of LIO for more de- 
tails). 

We used the Monte Carlo approach described in Broos et al. 
(2007, 2011) to estimate the false-match probability for each 
ONIR catalog. The main-catalog X-ray sources are consid- 
ered to consist of two populations: an "associated population" 
for which true counterparts are expected in an ONIR catalog, 
and an "isolated population" for which no counterparts are ex- 
pected (e.g., the true counterparts may be too faint or blended 
with other sources and thus not included) in an ONIR cata- 
log. We estimated the false-match probability for the associ- 
ated population by producing a mock ONIR counterpart for 
each X-ray source and running the likelihood-ratio matching 
procedure to find the counterpart recovery fraction. The off- 
set between the mock counterpart and the X-ray source is se- 
lected randomly based on the positional uncertainties, and the 
magnitude of the mock counterpart is drawn randomly from 
the expected magnitude distribution of the counterparts (de- 
rived previously in the likelihood-ratio matching procedure). 
The mock ONIR catalog is thus composed of the mock coun- 
terparts and the original ONIR catalog with source positions 
shifted and potential counterparts removed. To estimate the 
false-match probability for the isolated population, we shifted 
the X-ray source positions and recorrelated the shifted sources 
with the ONIR sources using likelihood-ratio matching. The 
above simulations were performed 100 times for each X-ray 
source population, and the results were used to solve for the fi- 
nal false-match probability for each ONIR catalog (see Broos 
et al. 2011 for details). The false-match probability for the 
associated population is generally smaller than that for the 
isolated population, and the final false-match probability for 
each ONIR catalog is < 4%. The expected mean false-match 
probability for the main-catalog sources is « 2.1%, derived 
by weighting the false-match probabilities of individual ONIR 
catalogs with the number of primary counterparts in each cat- 
alog. We note that the high identification rate, combined with 
the small false-match rate, provides independent evidence that 
the vast majority of our X-ray detections are robust. 

For the 24 main-catalog sources that do not have highly sig- 
nificant multiwavelength counterparts, we visually inspected 



the X-ray images and found that the majority of them have 
apparent or strong X-ray signatures. Of these 24 sources, 19 
were detected in the full band, with a median number of full- 
band counts of 49.8; 17 were detected in the soft band, with a 
median number of soft-band counts of 40.5; 9 were detected 
in the hard band, with a median number of hard-band counts 
of 58.7; and 17 were detected in at least two of the three stan- 
dard bands. We also investigated the Chandra events for these 
24 sources and concluded that they were not compromised by 
short-lived cosmic-ray afterglows. Of these 24 unidentified 
sources, 5 were previously detected in the L08 main catalog, 
3 were previously detected in the L08 supplementary CDF-S 
plus E-CDF-S Chandra catalog, and 16 were only detected in 
the 4 Ms observations. As for the nature of these 24 uniden- 
tified sources, we refer readers to § 4. 1 of LIO and references 
therein for detailed discussion of the possibilities. For exam- 
ple, 5 of these 24 unidentified sources are probably related 
to off-nuclear X-ray sources associated with nearby galaxies 
(e.g., Hornschemeier et al. 2004; Lehmer et al. 2006; note 
that, in this paper, we did not attempt a thorough identifica- 
tion of off-nuclear X-ray sources). 

4.4. Main-Catalog Details 

We summarize in Table 2 the columns (a total of 79) in the 

main Chandra X-ray source catalog; the main catalog itself is 
presented in Table 3. The details of the 79 columns are given 
below. 

1. Column 1 gives the source sequence number (i.e., XID). 
We list sources in order of increasing right ascension. 

2. Columns 2 and 3 give the right ascension and declination 
of the X-ray source, respectively. We determined source posi- 
tions following the procedure detailed in § 3.2. To avoid trun- 
cation error, we quote the positions to higher precision than in 
the International Astronomical Union (lAU) registered names 
that begin with the acronym "CXO CDFS". 

3. Columns 4 and 5 give the minimum value of \ogP (P 
is the AE-computed binomial no-source probability) among 
the three standard bands, and the logarithm of the mini- 
mum WAVDETECT false-positive probability detection thresh- 
old, respectively. More negative values of \ogP (Column 4) 
and false-positive probability threshold (Column 5) indicate 
a more significant source detection. We set logP = -99.0 for 
sources with P = 0. For the main-catalog sources, the me- 
dian value of logP is -8.9 (note that P < 0.004, correspond- 
ing to logP < -2.4, is the condition for a source to be in- 
cluded in the main catalog). There are 493, 57, 69, and 121 
sources with minimum WAVDETECT probabiHties^^ of 10"^, 
10"^, 10"^, and 10"^, respectively (see Fig. 5). 

4. Column 6 gives the Rs 68% confidence-level X-ray po- 
sitional uncertainty in arcseconds computed using eq. (2), 
which is dependent on both off-axis angle and aperture- 
corrected net source counts. The « 68% confidence-level 
X-ray positional uncertainty was used in the likelihood-ratio 
matching procedure (see § 4.3). The positional uncertainty for 
the main-catalog sources ranges from 0.10" to 1.51", with a 
median value of 0.42". 

5. Column 7 gives the off-axis angle of the X-ray source in 
arcminutes, which is the angular separation between the X-ray 
source (coordinates given in Columns 2 and 3) and the CDF-S 
average aim point (given in Table 1). The off-axis angle for 
the main-catalog sources ranges from 0.33' to 12.36', with a 
median value of 5.82'. The maximum off-axis angle of 12.36' 
is slightly larger than a half of the diagonal size of the ACIS-I 
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TABLE 2 

Overview of Columns in the Main Chandra Source Catalog 



Column Description 



1 Source sequence number (i.e., XID) 

2, 3 Right ascension and declination of the X-ray source 

4 Minimum value of log P among the three standard bands (P is the AE-computed binomial no-source probability) 

5 Logarithm of the minimum WAVDETECT false-positive probability detection threshold 

6 ~ 68% confidence-level X-ray positional uncertainty 

7 Off-axis angle of the X-ray source 

8-16 Aperture-corrected net (i.e., background-subtracted) source counts and the corresponding errors for the three standard bands 

17 Hag of whether a source shows any evidence for spatial extent 

18, 19 Right ascension and declination of the optical/near-infrared/infrared/radio (ONIR) counterpart 

20 Offset between the X-ray source and ONIR counterpart 

21 AB magnitude of the ONIR counterpart 

22 Name of the ONIR catalog from which the primary counterpart has been taken 

23^3 Right ascension, declination, and AB magnitude of the counterpart in seven ONIR catalogs 

A4-A6 Spectroscopic redshift, redshift quality flag, and the reference for the redshift 

47-57 Photometric-redshift information taken from sources in the literature 

58 Preferred redshift adopted in this paper 

59 Corresponding 2 Ms CDF-S source number from the main and supplementary Chandra catalogs presented in LOS 
60, 61 Right ascension and declination of the corresponding LOS source 

62 Corresponding 250 ks E-CDF-S source number from the main and supplementary Chandra catalogs presented in L05 

63, 64 Right ascension and declination of the corresponding LOS E-CDF-S source 

65-67 Effective exposure times derived from the exposure maps for the three standard bands 

68-70 Band ratio and the coiresponding errors 

71-73 Effective photon index with the corresponding errors 

74-76 Observed-frame fluxes for the three standard bands 

77 Absorption-corrected, rest-frame 0.5-8 keV luminosity 

78 Estimate of likely source type 

79 Notes on the source 



TABLE 3 
Main Chandra Source Catalog 



X-ray Coordinates Detection Probability Counts 



No. 


0:2000 


<520OO 


logP 


WAVDETECT 


Pos Err 


Off-Axis 


FB 


FB Upp Err 


FB Low Err 


SB 


SB Upp Err 


SB Low Err 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


1 ... 


03 31 35.79 


-27 51 36.0 


-99.0 


-8 


0.5 


11.98 


186.8 


19.0 


17.9 


117.8 


13.5 


12.4 


2... 


03 31 40.12 


-27 47 46.6 


-30.9 


-8 


0.5 


10.62 


155.7 


19.S 


1S.6 


101.5 


13.4 


12.2 


3 ... 


03 31 41.01 


-27 44 34.7 


-15.6 


-8 


0.6 


11.10 


96.5 


15.7 


14.5 


31.5 


S.5 


7.3 


4... 


03 31 43.25 


-27 54 05.6 


-6.5 


-5 


0.8 


11.41 


54.1 


13.8 


12.6 


19.9 


-1.0 


-1.0 


5... 


03 31 43.42 


-27 51 03.8 


-5.9 


-8 


0.5 


10.21 


109.0 


27.5 


25.7 


38.1 


14.7 


12.9 



Note. — Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminutes, and arcseconds. Table 3 is presented in its entirety in the electronic 
edition. A portion is shown here for guidance regarding its form and content. The full table contains 79 columns of information for the 740 X-ray sources. 



field of view (11.95'), due to the fact that the CDF-S obser- 
vations have varying aim points and roll angles, as shown in 
Table 1. 

6. Columns 8-16 give the aperture-corrected net (i.e., 
background-subtracted) source counts and the corresponding 
1(7 upper and lower statistical errors (Gehrels 1986) for the 
three standard bands, respectively. The photometry was cal- 
culated by AE using the position given in Columns 2 and 3 
for all bands and following the procedure described in § 3.2, 
and was not corrected for vignetting or exposure time vari- 
ations. To be consistent with our source detection criterion 
(i.e., P < 0.004), we considered a source to be "detected" 
for photometry purposes in a given band only if the AE- 
computed binomial no-source probability for that band is less 
than 0.004. For sources not detected in a given band, we cal- 
culated upper limits and placed -1.00 in the corresponding 
error columns. When the total number of counts within the 
polygonal extraction region of an undetected source was < 10, 
we computed the upper limit using the Bayesian method of 
Kraft et al. (1991) for a 99% confidence level; otherwise, we 



computed the upper hmit at the 3cr level for Poisson statistics 
(Gehrels 1986). 

7. Column 17 gives a flag indicating whether a source 
shows any evidence for spatial extent in basic testing. In 
§ 3.2, we ran wavdetect using 9 wavelet scales up to 16 
pixels, which potentially allows detection of sources that are 
extended on such scales. We utilized the following procedure 
to assess extent. We first derived a set of cumulative EEFs by 
extracting the PSF power within a series of circular apertures 
(centered at the source position) up to a 90% EEF radius from 
the merged PSF image. We then derived another set of cu- 
mulative EEFs by extracting source counts within a series of 
circular apertures (also centered at the source position) up to 
the same 90% EEF radius from the merged source image. Fi- 
nally, we used a Kolmogorov-Smirnov (K-S) test suitable for 
two distributions to compute the probability (pks) that the two 
sets of cumulative EEFs are consistent with each other. Of the 
740 main-catalog sources, 7 have pks < 0.01 (i.e., the merged 
PSF and source images are inconsistent with each other at or 
above a 99% confidence level) and have the value of this col- 
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umn set to 2; 24 have 0.01 < pks < 0.05 and have the value of 
this column set to 1 ; all the remaining sources have the value 
of this column set to 0. A total of 3 1 main-catalog sources are 
flagged as 1 or 2 that corresponds to a > 95% confidence level, 
which is comparable to the expected number of false-positive 
determinations, i.e., 37 =740 x (1—95%). These 31 sources 
are located across the entire CDF-S field and do not show the 
hkely expected pattern of central clustering (since the PSF is 
sharpest near the field center), which might also indicate that 
many of these sources could be false positives. Moreover, we 
did not find any significant signature of extension for these 3 1 
sources upon visual inspection. For the sources that truly have 
slight extents or are point sources sitting on top of highly ex- 
tended sources, our AE-computed photometry should be rea- 
sonably accurate, as detailed in § 3.2. We note that a few 
highly extended sources in the CDF-S (e.g., Giacconi et al. 
2002; LOS) cannot be identified here because these sources 
have larger extents than the maximum value of our adopted 
wavelet scales (i.e., 16 pixels); a full study of such extended 
sources is beyond the scope of this paper and will be presented 
in A. Finoguenov et al. (in preparation). 

8. Columns 18 and 19 give the right ascension and dec- 
lination of the ONIR counterpart (see § 4.3 for the details of 
multiwavelength identifications). Sources without multiwave- 
length identifications have these right ascension and declina- 
tion values set to "0 00 00.00" and "-00 00 00.0". 

9. Column 20 gives the measured offset between the X-ray 
source and ONIR counterpart in arcseconds. Sources without 
multiwavelength identifications have a value set to -1 .00. 

10. Column 21 gives the AB magnitude of the ONIR 
counterpart, measured in the counterpart-detection band.^'^ 
Sources without counterparts have a value set to -1 .00. 

11. Column 22 gives the name of the ONIR catalog (i.e., 
VLA, GOODS-S, GEMS, MUSIC, WFl, MUSYC, or SIM- 
PLE) from which the primary counterpart has been taken. 
Sources without counterparts have this column set to "...". 

12. Columns 23^3 give the right ascension, declination, 
and AB magnitude of the counterpart in the above seven 
ONIR catalogs that are used for identifications (i.e., WFI, 
GOODS-S, GEMS, MUSIC, MUSYC, SIMPLE, and VLA). 
We cross matched the positions of primary ONIR counter- 
parts (i.e.. Columns 17 and 18) with the seven ONIR cata- 
logs using likelihood-ratio matching. Sources without coun- 
terparts have corresponding right ascension and declination 
values set to "0 00 00.00" and "-00 00 00.0" and AB magni- 
tudes set to -1.00. We find « 75%, 61%, 72%, 55%, 70%, 
88%, and 18% of the main-catalog X-ray sources have WFI, 
GOODS-S, GEMS, MUSIC, MUSYC, SIMPLE, and VLA 
counterparts,"*-' respectively, with a false-match probability of 
< 2% for each ONIR catalog (see § 4.3 for details). 

13. Colunms 44^6 give the spectroscopic redshift (zspec). 
redshift quality flag, and the reference for the redshift. Spec- 
troscopic redshifts were collected from Le Fevre et al. (2004), 
Szokoly et al. (2004), Zheng et al. (2004), Mignoli et al. 
(2005), Ravikumar et al. (2007), Vanzefla et al. (2008), 
Popesso et al. (2009), Treister et al. (2009)'*'*, Balestra 

The AB magnitudes for the radio counterparts were converted from the 
radio flux densities, m(AB) =-2.51og(/i/)-48.60. 

'^^ Note that the GOODS-S and MUSIC catalogs cover « 39% of the 
CDF-S while the other five catalogs cover the entire CDF-S (see Table 1 
of LIO for more details); ~ 70% of the main-catalog sources are in the 
GOODS-S/MUSIC area [see Fig. 13(a)]. 

We flagged the spectroscopic redshifts from Treister et al. (2009) as 
"Insecure" since Treister et al. (2009) did not provide redshift quality flags. 



et al. (2010), and Silverman et al. (2010) with the ref- 
erence numbers of 1-10 in Column 46, respectively. We 
cross matched the positions of primary ONIR counterparts 
(i.e.. Columns 18 and 19) with the above catalogs of spectro- 
scopic redshifts using a matching radius of 0.5". Of the 716 
main-catalog sources that have multiwavelength identifica- 
tions, 419 (58.5%) have spectroscopic redshift measurements. 
343 (81.9%) of these 419 spectroscopic redshifts are secure, 
i.e., they are measured at > 95% confidence levels with mul- 
tiple secure spectral features (flagged as "Secure" in Column 
45); 76 (18.1%) of these 419 spectroscopic redshifts are in- 
secure (flagged as "Insecure" in Column 45). We estimated 
the false-match probability to be < 1 % in all cases. Sources 
without spectroscopic redshifts have these three columns set 
to -1.000, "None", and -1, respectively. 

14. Columns 47-57 give the photometric-redshift (zphot) 
information taken from sources in the literature. Columns 
47-50 give the photometric redshift, the corresponding Icr 
lower and upper bounds, and the alternative photometric 
redshift (set to -1.000 if not available) from LIO. Columns 
51-54 give the photometric redshift, the corresponding la 
lower and upper bounds, and the corresponding quality flag 
(smaller values of indicate better quality; < < 1-3 
indicates a reliable photometric-redshift estimate) from Car- 
damone et al. (2010). Columns 55-57 give the photometric 
redshift and the corresponding 1 a lower and upper bounds^^ 
from Rafferty et al. (201 1). We chose the above photometric- 
redshift catalogs because they utilized extensive multiwave- 
length photometric data and produced accurate photometric 
redshifts. LIO derived high-quality photometric redshifts for 
the 462 L08 main-catalog X-ray sources with a treatment of 
photometry that included utiUzing likelihood-matching, man- 
ual source deblending, and appropriate upper limits. Car- 
damone et al. (2010) employed new medium-band Subaru 
photometry and a PSF-matching technique to create a uni- 
form photometric catalog and derived photometric redshifts 
for over 80,000 sources in the E-CDF-S; their photometric 
redshifts are of high quality, in particular for bright sources. 
Rafferty et al. (2011) derived photometric redshifts for over 
100,000 sources in the E-CDF-S, using a compiled photomet- 
ric catalog that probes fainter magnitudes than the Cardamone 
et al. (2010) catalog by including sources in the GOODS-S 
MUSIC catalog (Grazian et al. 2006; Santini et al. 2009); 
their photometric redshifts are accurate down to faint fluxes. 
We cross matched the positions of primary ONIR counterparts 
(i.e.. Columns 18 and 19) with the above photometric-redshift 
catalogs using a matching radius of 0.5". Of the 716 main- 
catalog sources that have multiwavelength identifications, 668 
(93.3%) have photometric-redshift estimates from at least one 
source (this number excludes sources identified as stars, given 
in Column 78, that have all these columns set to -1.000). 
We estimate the false-match probability to be < 1% in all 
cases. Sources without photometric redshifts have aU these 
columns set to -1.000. We show in Fig. 8 the histograms of 
(a) (Zphot-Zspec)/(l+Zspec) and (b) Zphot for the above fliree 
sources of photometric redshift. It seems clear that the photo- 

The photometric-redshift errors derived with the Zurich Extragalac- 
tic Bayesian Redshift Analyzer (ZEBRA; Feldmann et al. 2006) gener- 
ally underestimate the real errors by factors of Ri 3 and « 6 for the spec- 
troscopic and non-spectroscopic samples, respectively (see, e.g., § 3.4 of 
LIO). Therefore, multiplying the photometric-redshift errors presented here 
by these corresponding factors (i.e., 3 and 6 for the spectroscopic and 
non-spectroscopic samples, respectively) wiU roughly give realistic la errors. 
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metric redshifts from each of these three sources have high 
quality"*'' in terms of accuracy and outlier percentage [see 
Fig. 8(a)] and cover a similar range of z « 0-5 [see Fig. 8(b)]. 
We refer readers to the cited references for the respective de- 
tails of the photometric-redshift derivations, the advantages of 
the adopted methodologies, and the caveats when using these 
photometric redshifts. 

15. Column 58 gives the preferred redshift adopted in 
this paper. We chose redshifts, in order of preference, 
as follows: (1) secure spectroscopic redshifts; (2) insecure 
spectroscopic redshifts that are in agreement with at least 
one of the LIO, Cardamone et al. (2010), or Rafferty 
et al. (2011) photometric-redshift estimates [i.e., |(Zspec- 
Zphot)/(l +Zspec)| < 0.15, where Zspec/Zphot is the spectro- 
scopic/photometric redshift]; (3) the LIO photometric red- 
shifts; (4) the Cardamone et al. (2010) photometric redshifts; 
and (5) the Rafferty et al. (201 1) photometric redshifts. Of the 
716 main-catalog sources that have multiwavelength identifi- 
cations, 673 (94.0%) have spectroscopic or photometric red- 
shifts. 

16. Column 59 gives the corresponding 2 Ms CDF-S source 
number from the main and supplementary Chandra catalogs 
presented in LOS. We matched our X-ray source positions 
(i.e.. Columns 2 and 3) to L08 source positions (corrected 
for the systematic positional shifts described in § 3.1) using 
a 2.5" matching radius for sources with off-axis angle 6 < 6' 
and a 4.0" matching radius for sources with > 6'. The mis- 
match probability is w 1% using this approach. For the 740 
main-catalog sources, we find 

(a) 440 have matches to the 462 LOS main-catalog sources 
(the value of Column 59 is that from col. [1] of Table 2 
in LOS; see § 4.5 for more details); 

(b) 41 have matches to the 86 LOS supplementary CDF-S 
plus E-CDF-S Chandra catalog sources (the value of 
Column 59 is that from col. [1] of Table 5 in LOS with 
a prefix of "SP1_", e.g., SP1_1); 

(c) 22 have matches to the 30 LOS supplementary opti- 
cally bright Chandra catalog sources (the value of Col- 
umn 59 is that from col. [1] of Table 6 in LOS with a 
prefix of "SP2_", e.g., SP2_1); 

(d) 6 were outside of the 2 Ms CDF-S footprint of LOS (the 
value of Column 59 is set to -1); the detection of these 
sources is simply due to the new sky coverage (rather 
than the improved sensitivity) of the 4 Ms CDF-S; and 

(e) 231 have no match in any of the LOS main and supple- 
mentary Chandra catalogs; these sources were inside 
the 2 Ms CDF-S footprint but are only detected now 
due to the improved sensitivity of the 4 Ms observations 
(the value of Column 59 is set to 0). 

In summary, of the 740 main-catalog sources, 503 were de- 
tected previously in the 2 Ms CDF-S observations (the value 
of Column 59 is greater than 0) and 237 were detected only 

'^^ In Fig. 8(a), the photometric redshifts fiom both LIO and Rafferty et al. 
(201 1) appear to have smaller outlier percentages than those from Cardamone 
et al. (2010) because the spectral energy distribution templates were opti- 
mized using the spectroscopic-redshift information before template fitting in 
both LID and Rafferty et al. (2011). Blind-test results show that the actual 
outlier percentages from LIO and Rafferty et al. (2011) are comparable to 
those from Cardamone et al. (2010) (see, e.g., § 3.4 of LIO for the details of 
blind tests). 



in the 4 Ms observations (the value of Column 59 is either 
-1 or 0). Compared to the LOS main catalog, there are 300 
(i.e., 740-440 = 300) new main-catalog sources (see § 4.7 for 
more details of these 300 sources). 

17. Columns 60 and 61 give the right ascension and decli- 
nation of the corresponding LOS source (corrected for the sys- 
tematic positional shifts described in § 3.1) indicated in Col- 
umn 59. Sources without an LOS match have right ascension 
and decHnation values set to "0 00 00.00" and "-00 00 00.0". 

18. Column 62 gives the corresponding 250 ks E-CDF-S 
source number from the main and supplementary Chandra 
catalogs presented in L05. We adopted the same matching ap- 
proach between X-ray catalogs as used for Column 59, again 
with the E-CDF-S source positions corrected for the system- 
atic positional shifts described in § 3.1. For the 740 main- 
catalog sources, we find (1) 239 have matches in the E-CDF-S 
main Chandra catalog (the value of Column 62 is that from 
col. [1] of Table 2 in L05); (2) 5 have matches in the E-CDF-S 
supplementary optically bright Chandra catalog (the value of 
Column 62 is that from col. [ 1 ] of Table 6 in L05 with a prefix 
of "SP_", e.g., SP_1); and (3) 496 have no match in either of 
the E-CDF-S main or supplementary Chandra catalogs (the 
value of Column 62 is set to 0). 

19. Columns 63 and 64 give the right ascension and decli- 
nation of the corresponding L05 E-CDF-S source (corrected 
for the systematic positional shifts described in § 3.1) indi- 
cated in Column 62. Sources without an E-CDF-S match have 
right ascension and declination values set to "0 00 00.00" and 
"-00 00 00.0". 

20. Columns 65-67 give the effective exposure times de- 
rived from the exposure maps (detailed in § 3.1) for the full, 
soft, and hard bands. Dividing the counts in Columns 8- 
16 by the corresponding effective exposure times will pro- 
vide effective count rates that have been corrected for vi- 
gnetting, quantum-efficiency degradation, and exposure time 
variations. 

21. Columns 68-70 give the band ratio and the corre- 
sponding upper and lower errors, respectively. We defined 
the band ratio as the ratio of counts between the hard and soft 
bands, correcting for differential vignetting between the hard 
and soft bands using the appropriate exposure maps. We fol- 
lowed the numerical error-propagation method described in 
§1.7.3 of Lyons (1991) to compute band-ratio errors. This 
method avoids the failure of the standard approximate vari- 
ance formula when the number of counts is small and the 
error distribution is non-Gaussian (e.g., see §2.4.5 of Eadie 
et al. 1971). We calculated upper limits for sources detected 
in the soft band but not the hard band and lower limits for 
sources detected in the hard band but not the soft band. For 
these sources, we set the upper and lower errors to the com- 
puted band ratio. We set band ratios and corresponding errors 
to -1 .00 for sources detected only in the full band. 

22. Columns 71-73 give the effective photon index (F) with 
the corresponding upper and lower errors, respectively, for 
a power-law model with the Galactic column density given 
in § 1. We calculated the effective photon index based on 
the band ratio in Column 68, using a conversion between the 
effective photon index and the band ratio. We derived this 
conversion using the band ratios and photon indices calcu- 
lated by the AE-automated XSPEC-fitting procedure for rel- 
atively bright X-ray sources (with full-band counts greater 
than 200; this ensures rehable XSPEC-fitting results). This 
approach takes into account the multi-epoch Chandra caU- 
bration information and thus has an advantage over meth- 
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Fig. 8. — (a) Histogram of (Zphot ~Zspec)/(l +Zspec) for LIO (218 sources), Cardamone et al. (2010; 314 sources), and Rafferty et al. (201 1; 339 sources), (b) 
Histogram of Zphot for LIO (black histogram; 417 sources), Cardamone et al. (2010; dark-gray histogram; 508 sources), and Rafferty et al. (2011; hght-gray 
histogram; 611 sources). The histograms have been slightly shifted for clarity. 



ods using only single-epoch calibration information such as 
the CXC's Portable, Interactive, Multi-Mission Simulator 
(PIMMS) method used by LOS. We calculated upper limits 
for sources detected in the hard band but not the soft band 
and lower limits for sources detected in the soft band but not 
the hard band. For these sources, we set the upper and lower 
errors to the computed effective photon index. For low-count 
sources, we are unable to determine the effective photon index 
reliably; we therefore assumed F = 1 .4, which is a representa- 
tive value for faint sources that should yield reasonable fluxes, 
and set the corresponding upper and lower errors to 0.00. We 
defined sources with a low number of counts as being (1) de- 
tected in the soft band with < 30 counts and not detected in 
the hard band, (2) detected in the hard band with < 15 counts 
and not detected in the soft band, (3) detected in both the soft 
and hard bands, but with < 15 counts in each, or (4) detected 
only in the full band. 

23. Columns 74-76 give observed-frame fluxes in units of 
erg cm~^ s"' in the full, soft, and hard bands. We computed 
fluxes using the counts in Columns 8-16, the appropriate ex- 
posure maps (Columns 65-67), and the effective power-law 
photon indices given in Column 7 1 . We did not correct fluxes 
for absorptions by Galactic material or material intrinsic to the 
source. Negative flux values indicate upper limits. We note 
that, due to the Eddington bias, sources with low net counts 
(given in Columns 8-16) could have true fluxes lower than 
those computed here (see, e.g., VikhHnin et al. 1995; Geor- 
gakakis et al. 2008). We do not attempt to correct for the 
Eddington bias, since we aim to provide only observed fluxes 
here. Determining more accurate fluxes for these sources 
would require (1) using a number-count distribution prior to 
estimate the flux probabilities for sources near the sensitivity 
limit and/or (2) directly fitting the X-ray spectra for each ob- 
servation; these analyses are beyond the scope of this paper. 

24. Column 77 gives a basic estimate of the absorption- 
corrected, rest-frame 0.5-8 keV luminosity (Lo.5-8 kev) in 
units of erg s"'. We calculated Lo.5-8 keV using the procedure 
detailed in § 3.4 of Xue et al. (2010). Briefiy, this procedure 
models the X-ray emission using a power-law with both in- 
trinsic and Galactic absorption (i.e., zpow x wabs x zwabs in 
XSPEC) to find the intrinsic column density that reproduces 
the observed band ratio (given in Coluimi 68), assuming a typ- 



ical power-law photon index of Fim =1.8 for intrinsic AGN 
spectra. It then corrects for both Galactic and intrinsic ab- 
sorption to obtain the absorption-corrected flux (/o.s-g kev.int; 
as opposed to the observed flux given in Column 74), and fol- 
lows the equation Lo.5-8 kev = 47rc/^/o.5-8 kev,int(l +z)'"'°'"^ to 
derive Lo.5-8 keV (where di is the luminosity distance and z is 
the adopted redshift given in Column 58). In this procedure, 
we set the observed band ratio to a value that corresponds to 
r = 1 .4 for sources detected only in the full band; for sources 
having upper or lower Hmits on the band ratio, we adopted 
their upper or lower limits for this calculation. Basic lumi- 
nosity estimates derived in this manner are generally found to 
agree with those from direct spectral fitting to within a factor 
of ~ 30%"*^; the direct spectral-fitting approach should pro- 
duce more reliable estimates, but is beyond the scope of this 
paper. Sources without redshift estimates have this colunon 
set to -1.000; negative luminosity values other than -1.000 
indicate upper limits. 

25. Column 78 gives a basic estimate of likely source 
type. We categorized the X-ray sources into three basic types: 
"AGN", "Galaxy", and "Star". We utiUzed four criteria that 
are based on distinct AGN physical properties and one cri- 
terion that is based on optical spectroscopic information to 
identify AGN candidates, which must satisfy at least one of 
these five criteria. We briefly describe these criteria below: 

(a) A source with an intrinsic X-ray luminosity (given in 
Column 77) of Lo.5-8 keV > 3 x lO"*^ erg s"' wiU be iden- 
tified as a luminous AGN. 

(b) A source with an effective photon index (given in Col- 
umn 71) of r < 1.0 will be identified as an obscured 
AGN. 

(c) A source with an X-ray-to-optical flux ratio of 
log(/x//fi) > -1 (where fx = /0.5-8 kev, /0.5-2 kev, or 
fi-s kev) wiU be identified as an AGN. 

We caution that our basic io.5-8 kev estimates could be subject to larger 
uncertainties for heavily obscured AGNs. This is not only due to the increas- 
ing difficulty in determining the intrinsic column density from the observed 
band ratio, but also due to the fact that other components (e.g., reflection, and 
scattering) become stronger in such heavily obscured sources. 
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(d) A source with excess (i.e., a factor of > 3) X-ray emis- 
sion over the level expected from pure star formation 
will be identified as an AGN, i.e., with Lo.s-gkeV ^ 
3 X (8.9 X 10"Lr), where Lr is the rest-frame 1.4 GHz 
monochromatic luminosity in units of W Hz~' and 
8.9 X IO'^Lr is the expected X-ray emission level that 
originates from starburst galaxies (see Alexander et al. 
2005 for the details of this criterion). 

(e) A source with optical spectroscopic AGN features such 
as broad emission lines and/or high-excitation emis- 
sion lines will be identified as an AGN; we cross 
matched the sources (using the ONIR counterpart po- 
sitions given in Columns 18 and 19) with the spectro- 
scopically identified AGNs in Szokoly et al. (2004), 
Mignoli et al. (2005), and Silverman et al. (2010), us- 
ing a matching radius of 0.5". 

We note that the above five criteria are effective but not com- 
plete in identifying AGNs and refer readers to, e.g., Bauer 
et al. (2004), Alexander et al. (2005), Lehmer et al. (2008), 
and Xue et al. (2010) for discussions and caveats (e.g., low- 
luminosity and/or highly obscured AGNs may still not be 
identified through the criteria presented here). We also iden- 
tified likely stars by cross matching the sources (using the 
ONIR counterpart positions given in Columns 18 and 19) 
with (1) the spectroscopically identified stars in Szokoly et al. 
(2004), Mignoh et al. (2005), and Silverman et al. (2010), 
(2) the likely stars with stellarity indices greater than 0.7 
in the GEMS HST catalog (Caldwell et al. 2008), and (3) 
the likely stars with best-fit stellar templates in the MUSYC 
photometric -redshift catalog (Cardamone et al. 2010), using 
a matching radius of 0.5". We inspected each of the sources 
identified as stars in the 7/5r images and retrieved sources that 
appear to be galaxies (i.e., set our classification to galaxy). 
The sources that were not identified as AGNs or stars are clas- 
sified as "galaxies". Of the 740 main-catalog sources, 568 
(76.8%), 162 (21.9%), and 10 (1.3%) are identified as AGNs, 
galaxies, and stars, respectively. Of the 568 AGNs in the main 
catalog, 65.1%, 40.3%, 91.7%, 14.8%, and 1.1% satisfy the 
criteria (a), (b), (c), (d), and (e), respectively. 

26. Column 79 gives notes on the sources. We annotated 
sources at the field edge that lie partially outside of the sur- 
vey area with "E" (one source only) and sources in close 
doubles or triples with "C" (a total of 35 sources; these 35 
sources have overlapping polygonal extraction regions that 
correspond to « 40-75% EEFs; see § 3.2). Sources not an- 
notated have this column set to 

4.5. Comparison with 2 Ms CDF-S Main-Catalog Sources 

We summarize in Table 4 the source detections in the three 
standard bands. In total 740 sources are detected, with 634, 
650, and 403 detected in the full, soft, and hard band, respec- 
tively. As stated earlier in § 4.4 (see the description of Col- 
umn 59), 503 of the main-catalog sources were detected in 
the L08 main or supplementary catalogs, among which 440 
were detected in the L08 main catalog. For these 440 com- 
mon sources, we find general agreement between the derived 
X-ray photometry presented here and in L08. For instance, the 
median ratio between our full-band count rates and the L08 
full-band count rates for the 387 full-band detected sources 
(among these 440 common sources) is 1.03, with an interquar- 
tile range of 0.91-1.14. The « 3% increase in the full-band 
count rates is mainly caused by a few updates to the ancillary 



TABLE 4 

Summary of Chandra Source Detections 



Number of Detected Counts Per Source 

Band (keV) Sources Maximum Mmimum Median Mean 

Full (0.5-8.0) 634 35657.0 11.4 101.4 497.4 

Soft (0.5-2.0) 650 25470.7 6.0 45.1 293.8 

Hard (2-8) 403 10219.3 10.7 99.9 302.8 



response file (ARE) and contamination model in the CALDB 
data^** since the production of the L08 catalogs. The detailed 
differences (e.g., scattering) in the derived X-ray photometry 
are mainly due to source variability and/or the above caldb 
updates (e.g., sources with different X-ray spectral shapes are 
affected differently by these caldb updates). The approx- 
imately doubled exposure improves the source positions and 
spectral constraints significantly. Hence, the 4 Ms CDF-S cat- 
alogs presented here supersede those in L08. 

Twenty-two (i.e., 462-440 = 22) of the 462 sources de- 
tected in the L08 main catalog are not included in our main 
catalog, among which 3 are included in our supplementary 
catalog (see § 5). Thus, there are a total of 19 "missing" L08 
main-catalog sources not included in the 4 Ms main or sup- 
plementary catalogs. Among these 19 missing sources, there 
are two cases where a source was previously listed as being 
in a close pair but is now removed due to no apparent sig- 
nature of a close pair in the 4 Ms images. Of the remaining 
17 sources, 12, 3, and 2 have a logarithm of the minimum 
WAVDETECT false-positive probability detection threshold of 
-6, -7, and -8 in the L08 main catalog, respectively. Among 
these 17 sources, 9 have no multi wave length counterparts and 
have no emission clearly distinct from the background in the 
4 Ms images, which indicates that most of these 9 sources are 
likely false detections.^^ Eor the other 8 sources that have rel- 
atively faint multiwavelength counterparts, they also have no 
apparent X-ray signatures in the 4 Ms images although a few 
of them have full-band counts of > 20-30 in the L08 catalog; 
these 8 sources are likely real X-ray sources, but they are not 
detected in the 4 Ms images probably due to source variabil- 
ity and/or background fluctuations, as the second 2 Ms expo- 
sure was taken w 2.5 years after the completion of the first 
2 Ms exposure. Indeed, all of these 8 sources were variable 
(at > 99.7% confidence levels based on K-S tests) and became 
fainter (i.e., having a factor of > 2 smaller count rates) during 
the second 2 Ms of observations; consequently, the addition of 
background counts diluted their signals from the first 2 Ms of 
observations.^*^ We note that source variability is not uncom- 
mon among the CDF sources: over short timescales (days-to- 
weeks) a w 35% median flux variabiHty for the sources in the 
first 2 Ms data set has been observed; over long timescales 
(years) source fluxes could vary by up to a factor of w5-10 
in a few extreme cases (Paolillo et al. 2004; M. Paolillo et al. 
201 1, in preparation). 

We sunmiarize in Table 5 the number of sources detected 

For example, there was a recalibration of the ACIS-I ARE in 
CALDB 4. LI (released in January 2009), yielding a flat Ri 9% reduction in 
the effective area below 2 keV and a « 0-8% reduction between 2 and 5 keV 
(see http://cxc.harvard.edu/ciao/why/caldb4. 1 . l_hrma.html). 

''^ We note that L08 estimated the number of false detections in their main 
catalog to be 18, which is a conservative estimate; the real number of false 
detections is likely 2-3 times smaller, i.e., 6-9 (see § 3.2 of LOS). 

Of these 8 sources, only 3 satisfy the P < 0.004 criterion during the first 
2 Ms exposure, while none satisfies the P < 0.004 criterion during the second 
2 Ms exposure. 
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TABLE 5 

Sources Detected in One Band but not 
Another 



Detection Band Nondetection Energy Band 
(keV) F5II SoH HSd" 



Full (0.5-8.0) ... 85 236 

Soft (0.5-2.0) 101 ... 316 

Hard (2-8) 5 69 



Note. — For example, there were 85 
sources detected in the full band but not in 
the soft band. 
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Fig. 9. — Histograms of detected source counts for the main-catalog 
sources in the full (top), soft (middle), and hard (bottom) bands. Sources 
with upper limits have not been included in the plots. The vertical dotted 
lines indicate median numbers of counts in each band (see Table 4). 

in one band but not another. There are 21, 101, and 5 sources 
detected only in the full, soft, and hard band, respectively, as 
opposed to 3 1, 56, and 3 sources detected only in the full, soft, 
and hard band in the L08 main catalog. 

4.6. Properties of Main-Catalog Sources 

In Figure 9 we show the distributions of detected counts in 
the three standard bands for the main-catalog sources. The 
median number of counts is « 101, 45, and 100 for the full, 
soft, and hard band, respectively. There are 319 sources with 
> 100 full-band counts, for which basic spectral analyses are 
possible; there are 202, 101, and 60 sources with > 200, 500, 
1000 full-band counts, respectively. 

In Figure 10 we show the distributions of X-ray flux in the 
three standard bands for the main-catalog sources. The X-ray 
fluxes span roughly four orders of magnitude, with a median 
value of 6.8 x 10-'^ 1.0 x IQ-^^, and 1.1 x lO''^ erg cm'^ s'^ 
for the full, soft, and hard band, respectively. 

We show in Figure 1 1 the distribution of the AE-computed 
no-source probability P (given in Column 4) for the main- 
catalog sources; sources without multiwavelength counter- 
parts (given in Columns 18 and 19) are highlighted by shaded 
areas. It is clear that the majority of the main-catalog sources 
have low no-source probabilities (i.e., with logP < -6). We 
find that 1.3% of the logP < -6 sources have no multiwave- 
length counterparts, as opposed to the 6.6% of logP > -6 
sources that lack multiwavelength counterparts. Combined 
with the smaU false-match rate (see § 4.3), the above observa- 
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Fig. 10. — Histograms of X-ray fluxes for the main-catalog sources in the 
full (top), soft (middle), and hard (bottom) bands. Sources with upper limits 
have not been included in the plots. The vertical dotted lines indicate the 
median fluxes of 6.8 X 10"'*, 1.0 X 10"'* and 1.1 X 10"'^ erg cm"^ s"^ for 
the full, soft, and hard bands, respectively. 
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Fig. 1 1 . — Histogram of the AE-computed binomial no-source probability, 
P, for the main-catalog sources. For the purpose of illustration, we set the 
values of P < 10"^" to P = lO^^" in this plot. The shaded areas indicate 
sources that have no multiwavelength counterparts, with the numbers of these 
imidentified sources listed above the corresponding shaded areas. 

tions suggest that an X-ray source having a secure multiwave- 
length counterpart is an effective indicator of it being real. 

We show in Figure 12 "postage-stamp" images from the 
WFl R-hund, the GOODS-S/GEMS HST z-band, and the 
SIMPLE IRAC 3.6 //m band with adaptively smoothed fuU- 
band X-ray contours overlaid for the main-catalog sources. 
The size of X-ray sources in these images spans a wide range 
largely due to PSF broadening with off-axis angle. 

4.7. Properties of the 300 New Main-Catalog Sources 

In this section we examine the properties of the 300 main- 
catalog sources that were not detected in the L08 main catalog 
(hereafter new sources), putting emphasis on the comparison 
with the sources previously detected in the L08 main catalog 
(hereafter old sources; a total of 740 - 300 = 440 sources). 

Figure 13(a) shows the positions of the new sources (shown 
as filled symbols) and the old sources (shown as open sym- 
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Fig. 12.— Typical postage-stamp images from (a) the WFI R-band, (b) the GOODS-S/GEMS HST z-band, and (c) the SIMPLE IRAC 3.6 /im band for the 
main-catalog sources with full-band adaptively smoothed X-ray contours overlaid. The contours have a logarithmic scale and range from RiO.003%-30% of the 
maximum pixel value. The labels at the top of each image give the source name (for right ascension, the hours "03" have been omitted for succinctness) derived 
from the source coordinates and the source type ("A" denotes "AGN"; "G" denotes "Galaxy"; and "S" denotes "Star"). The numbers at the bottom of each image 
indicate the source number, the adopted redshift, and the full-band counts or upper limit (with a "<" sign). There are several cases where no X-ray contours are 
present, either because these sources were not detected in the full band or their full-band counts are low resulting in their observable emission in the adaptively 
smoothed images being suppressed by csmooth. Each image is 25" x 25", with the source of interest located at the center. The cutouts of all the main-catalog 
sources are available in the electronic edition. \See the electronic edition of the Supplement for a complete version of this figure.] 
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bols), with source types (given in Column 78) being color- 
coded (red for AGNs, black for galaxies, and blue for stars, 
respectively). Different symbol sizes represent different AE 
binomial no-source probabilities (see Column 4 of Table 3), 
with larger sizes indicating lower no-source probabilities (i.e., 
higher source-detection significances). In the GOODS-S re- 
gion, there are 512 main-catalog sources, with 221 being new; 
in the CANDELS region, there are 258 main-catalog sources, 
with 123 being new; and in the UDF region, there are 45 
main-catalog sources, with 20 being new. The source densi- 
ties of both new and old sources decline toward large off-axis 
angles as the sensitivity decreases (see § 7.2); such a trend 
appears more apparent among new sources than among old 
sources, e.g., 22.0% of new sources and while only 14.8% of 
old sources have < 3', and 62.0% of new sources and while 
only 46.1% of old sources have 9 < 6', respectively. Fig- 
ure 13(c) presents the observed source density as a function of 
off-axis angle for all the main-catalog sources. Overall, AGNs 
have larger observed source densities than galaxies. However, 
since the slope of the observed galaxy number counts at faint 
fluxes is steeper than that of the observed AGN number counts 
(e.g., Bauer et al. 2004), the galaxy source density approaches 
the AGN source density toward smaller off-axis angles (i.e., 
toward lower flux levels). This can also be seen in Fig. 13(d) 
that plots the observed source density versus off-axis angle 
for new sources; within 9 = 3', the new galaxies already out- 
number the new AGNs (36 versus 30). Near the center of the 
4 Ms CDF-S (within 6 = 3'), as shown in Fig. 13(c), the over- 
all observed AGN and galaxy source densities have reached 
9800!|^!!|| deg-2 and 6900!',',?' deg^^ respectively. We note 
that detailed analyses of the overall source densities for differ- 
ent source types, which consider effects such as the Eddington 
bias and incompleteness, are beyond the scope of this work. 

We show in Figure 14 plots of (a) observed-frame full- 
band flux (given in Column 74) vs. redshift (given in Col- 
umn 58), (b) absorption-corrected, rest-frame 0.5-8 keV lu- 
minosity (given in Column 77) vs. redshift, and (c) band ra- 
tio (given in Column 68) vs. absorption-corrected, rest-frame 
0.5-8 keV luminosity, for new sources (shown as filled cir- 
cles) and old sources (shown as open circles), respectively. 
Compared to old sources, new sources typically have smaller 
full-band fluxes and 0.5-8 keV luminosities [see their cluster- 
ing at the faint-flux end in Fig. 14(a) and at the low-luminosity 
end in Fig. 14(b)], which is expected since the 4 Ms CDF-S 
has fainter flux limits than the 2 Ms CDF-S. The existence 
of a small number of new sources at the high-flux/luminosity 
end leads to the full range of flux/luminosity for new sources 
being similar to that for old sources; these bright/luminous 
sources are typically located at relatively large off-axis an- 
gles. As shown in Fig. 14(a), there is no apparent correla- 
tion between full-band flux and redshift for either new or old 
sources, and the 4 Ms CDF-S is detecting an appreciable num- 
ber of the faintest sources at least up to z « 3. According 
to Fig. 14(b), the 0.5-8 keV luminosity spans a very broad 
range (roughly six orders of magnitude) for both new and old 
sources; 13.6% of the main-catalog sources are very lumi- 
nous (with Lo.5-8 kev > 10^ erg s"'; most are old sources), 
among which there are a number of sources that are highly 
obscured [see the upper right corner of Fig. 14(c)]. As seen 
in Fig. 14(c), new sources could potentially have a similar 
range or distribution of band ratio to that of old sources, given 
that 82.7% of new sources have either lower limits (19.0%) or 
upper limits (8 1 .0%) on their band ratios (see relevant discus- 



sions on this point later in this section). 

Figure 15 shows distributions of observed-frame full -band 
flux (given in Colurmi 74) and absorption-corrected, rest- 
frame 0.5-8 keV luminosity (given in Column 77) for new 
sources (main panels) and old sources (insets), separated by 
source type. Based on our source-classification scheme, it 
is clear that sources with different types have disparate dis- 
tributions of flux and luminosity when either new or old 
sources are considered, and that overall galaxies become the 
numerically dominant population at fuU-band fluxes less than 
« 10~'^ erg cm"-^ s~' or 0.5-8 keV luminosities less than 
« 10^^ erg s"';^' this trend is more pronounced when only 
new sources are considered. It is also clear that (1) new 
sources (either AGNs or galaxies) have similar ranges of flux 
and luminosity to those of old sources (either AGNs or galax- 
ies); (2) new sources (either AGNs or galaxies), as expected, 
typically have smaller fluxes (i.e., have a smaller median flux) 
than old sources (either AGNs or galaxies); (3) compared 
to old AGNs, new AGNs typically have smaller luminosities 
(i.e., have a smaller median luminosity); and (4) compared to 
old galaxies, new galaxies have comparable luminosities (i.e., 
have about the same median luminosity). 

We show in Figure 16(a) the band ratio as a function of full- 
band count rate for new sources (shown as filled symbols) 
and old sources (shown as open symbols). The sources are 
color-coded according to their likely types, with red, black, 
and blue colors indicating AGNs, galaxies, and stars, respec- 
tively. Also shown in Figure 16(a) are the average band ratios 
derived from stacking analyses following the procedure de- 
scribed in Luo et al. (201 1), for all AGNs, all galaxies, and all 
sources (including both AGNs and galaxies), shown as large 
crosses, triangles, and diamonds, respectively. As expected, 
the overall average band ratio is dominated by AGNs because 
most of the main-catalog sources are AGNs and AGNs typ- 
ically are more X-ray luminous than galaxies (see Fig. 15). 
The overall average band ratio rises between full-band count 
rates of « 10"^ and « 10""^ count s"\ and it levels off and 
subsequently decreases below full-band count rates of ~ 10"^ 
count s"' . The former increasing trend of the average band ra- 
tio is due to an increase in the number of absorbed AGNs de- 
tected at fainter fluxes and has been reported previously (e.g., 
Tozzi et et. 2001; A03; L05; L08); the latter decreasing trend 
of the average band ratio is partly because the contribution 
from normal and starburst galaxies increases at these lowest 
count rates (e.g., Bauer et al. 2004). Note that, at the low- 
est count rates studied, most of the sources have only band- 
ratio upper limits; thus the average band ratio lies below the 
individual-source upper limits. We show in Figure 16(b) the 
fraction of new sources as a function of fuU-band count rate 
for the main-catalog sources. Above full-band count rates 
of w 10"^ count s"', the fraction of new sources is small 
and roughly constant (~ 5-13%); below full-band count rates 
of « 10"^ count s~\ the fraction of new sources rises from 
« 12% to w 67% toward smaller full-band count rates. 

To examine further the band-ratio behavior of new and old 
sources, we show in Fig. 17 the average (i.e., stacked) band 

^' There may be a selection effect that can potentially contribute to 
the result that galaxies numerically dominate over AGNs at Lo.5-8 keV < 

10*- erg s"' since we used Lo.5-8 keV > 3 X 10''^ erg s"' as one of the AGN 
identification criteria. However, as shown in § 4.4 (see the description of 
Column 78), > 92% of the AGNs in the main catalog can be identified by the 
criteria other than the luminosity criterion; therefore, such a selection effect 
should be minimal. 
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Fig. 13. — (Top) Source spatial distributions for (a) the main catalog and (b) the supplementary optically bright catalog. Sources that are considered AGNs, 
galaxies, and stars (given in Column 78) are colored red, black, and blue, respectively. Open circles indicate AGNs/galaxies that were previously detected in (a) 
the LOS main catalog or (b) the LOS main or supplementary optically bright catalog; open stars in (a) indicate stars that were previously detected in the LOS main 
catalog; filled circles and stars indicate new AGNs/galaxies and stars, respectively. The regions and the plus sign are the same as those in Fig. 1 . In panel (a), 
the sizes of the circles and stars indicate the AE binomial no-source probabilities, with larger sizes indicating lower no-source probabilities: as the size becomes 
smaller, the AE binomial no-source probability P moves from logP < -5, -5 < logP < -4, -4 < logP < -3, to logP > —3. In panel (b), all sources have 
logP > -3 and are plotted as circles/stars of the same size. (Bottom) Observed source density for different source types as a function of off-axis angle for (c) 
all the main-catalog sources and (d) the new main-catalog sources, as computed in bins of A6 = 1'. la errors are calculated using Poisson statistics, [see the 
electronic edition of the Supplement for a color version of this figure.] 
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Fig. 14. — Plots of (a) observed-frame full-band flux vs. redshift, (b) absoiption-corrected, rest-frame 0.5-S keV luminosity vs. redshift, and (c) band ratio vs. 
absorption-corrected, rest-frame 0.5-8 keV luminosity, for the main-catalog sources. Gray open circles indicate the main-catalog sources that were previously 
detected in the LOS main catalog; black filled circles indicate the new main-catalog sources that were not previously detected in the LOS main catalog. Arrows 
indicate limits. Several sources shown in Panels (a) and (b) have photometric redshifts greater than ft: 4.5; these photometric redshifts are probably not veiy 
rehable due to poor photometric coverage (see § 3.3 of LIO for more discussion). In Panel (b), sources without redshift estimates have not been included in the 
plot; in Panel (c), sources without redshift estimates and sources with only full-band detections have not been included in the plot. The dotted lines in Panels (b, 
c) and the dashed-dot line in Panel (c) indicate the threshold values of two AGN-identification criteria (i.e., Z.o.5-8 kcV 3 X 10"*- erg s"' and F < 1.0; see the 
description of Column 7S for details). 
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Fig. 15. — Distributions of (a) observed-frame full-band flux and (b) 
absorption-corrected, rest-frame 0.5-8 keV luminosity for the new main- 
catalog sources. The red and black histograms indicate AGNs and galax- 
ies, respectively. The vertical red and black dashed lines indicate the median 
values for AGNs and galaxies, respectively. Sources with upper limits on 
full-band fluxes have not been included in the plotting for panel (a); sources 
without estimates of X-ray luminosities (due to no available redshift) or with 
upper limits on X-ray luminosities have not been included in the plotting for 
panel (b). The insets show results for the old main-catalog sources, [see the 
electronic edition of the Supplement for a color version of this figure.} 

ratio in bins of redshift and X-ray luminosity for new AGNs, 
old AGNs, new galaxies, and old galaxies, respectively. Ac- 
cording to Fig. 17, (1) new AGNs have larger band ratios than 
old AGNs no matter which bin of redshift or X-ray luminosity 
is considered (presumably due to the fact that the detection of 
highly absorbed AGNs with large band ratios requires deep 
observations given the small ACIS-I effective area at high 
energies); (2) in the two lower redshift bins (0 < z < 1 and 
1 < z < 2), the band ratios of new and old galaxies appear 
roughly consistent and constant within errors (hinting at no 
evolution in the X-ray spectral shape of the accreting binary 
populations that dominate the X-ray emission of normal and 
starburst galaxies); (3) new AGNs and old AGNs have similar 
patterns of band ratio versus X-ray luminosity, both peaking 
at the bin of 42.5 < log(Lx) < 43.5; and (4) in the lowest 
luminosity bin [log(Lx) < 41.5], new galaxies have a larger 
average band ratio than old galaxies, while in a higher lumi- 
nosity bin [41 .5 < log(Lx) < 42.5], new and old galaxies have 
consistent band ratios. 

We show in Figure 18(a) the WFI 7? -band magnitude ver- 
sus the full-band flux for new sources (filled symbols) and 
old sources (open symbols), as well as the approximate flux 
ratios for AGNs and galaxies (e.g., Maccacaro et al. 1988; 
Stocke et al. 1991; Homschemeier et al. 2001; Bauer et al. 



2004; also see the description of Column 78 for AGN iden- 
tification). The sources are color-coded according to their 
likely types, with red, black, and blue colors indicating AGNs, 
galaxies, and stars, respectively. For comparison, we also 
show in Fig. 18(c) the IRAC 3.6 fim magnitude versus the 
full-band flux for new sources (filled symbols) and old sources 
(open symbols), since a higher fraction of the main-catalog 
sources have counterparts in the IRAC 3.6 /im band than in 
the WFI R-hanA (i.e., « 88% vs. w 75%; see the description 
of Columns 23^3). Overall, a total of 568 (76.8%) of the 
main-catalog sources are likely AGNs, and the majority of 
them lie in the region expected for relatively luminous AGNs 
[i.e., log(/x//R) > -1; dark gray areas in Fig. 18(a)]; of these 
568 AGNs, 192 (33.8%) are new. A total of 162 (21.9%) of 
the main-catalog sources are likely galaxies, and the majority 
of them lie in the region expected for normal galaxies, star- 
burst galaxies, and low-luminosity AGNs [i.e., log(/x//R) < 
-1; light gray areas in Fig. 18(a)]; of these 162 sources, 104 
(64.2%) are new. Only 10 (1.3%) of the main-catalog sources 
are likely stars with low X-ray- to-optical flux ratios; of these 
10 stars, 4 are new. Among new sources, normal and starburst 
galaxies account for a fraction of w 35%, in contrast to w 13% 
if old sources are considered. The above source-classification 
analysis indicates that, as expected, the 4 Ms CDF-S survey 
is detecting sources in or close to a regime where the galaxy 
number counts approach the AGN number counts, due to the 
steeper number-count slope expected for galaxies (e.g., Bauer 
et al. 2004). 

Figure 19 shows the distributions of X-ray- to-optical flux 
ratio for new AGNs, old AGNs, new galaxies, and old galax- 
ies, respectively. It is clear that (1) new AGNs generally have 
smaller X-ray- to-optical flux ratios than old AGNs and (2) 
new and old galaxies have similar distributions of X-ray-to- 
optical flux ratio. 

5. SUPPLEMENTARY OPTICALLY BRIGHT CHANDRA SOURCE 
CATALOG 

5.1. Supplementary Catalog Production 

Of the 152 (i.e., 892-740 = 152) candidate-list X-ray 
sources that were not included in the main Chandra cata- 
log (see § 4.1), 119 were of moderate significance (0.004 < 
P < 0.1). To recover genuine X-ray sources from this sam- 
ple, we constructed a supplementary Chandra source catalog 
consisting of the subset of these sources that have bright op- 
tical counterparts. Given that the density of optically bright 
sources on the sky is comparatively low, it is likely that the 
Chandra sources with such counterparts are real. We selected 
bright optical sources from the WFI /?-band source catalog de- 
scribed in § 3.1. We searched for bright optical counterparts 
(having R < 24.0) to the 1 19 X-ray sources of interest using 
a matching radius of 1.2". The choices of 0.004 < P < 0.1, 
the /?-band cutoff magnitude, and the matching radius were 
made to ensure a good balance between the number of de- 
tected sources and the expected number of false sources. We 
find a total of 36 optically bright X-ray sources, of which 3 are 
L08 main-catalog sources that were not included in our main 
catalog and 3 are L08 supplementary optically bright sources 
(i.e., 30 new sources in the 4 Ms supplementary catalog). We 
note that the majority (22 out of 30) of the L08 supplemen- 
tary optically bright sources are included in our main catalog 
(see the description of Column 59); this explains the small 
number of L08 supplementary sources included in our sup- 
plementary catalog. We estimated the expected number of 
false matches to be w 2.2 (i.e., w 6.1%) by manually shifting 
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Fig. 16. — (a) Band ratio vs. full-band count rate for the main-catalog sources (for reference, the top ,v-axis shows representative full-band fluxes, which 
are derived from full-band count rates assuming a F = 1.4 power law). Sources that are considered AGNs, galaxies, and stars are colored red, black, and blue, 
respectively. Open circles and stars indicate AGNs/galaxies and stars that were previously detected in the LOS main catalog; filled circles and stars indicate 
new AGNs/galaxies and stars, respectively. Arrows indicate upper or lower limits, which mostly lie in the area of low count rates. Sources detected only 
in the full band cannot be plotted. Large crosses, triangles, and diamonds show average band ratios as a function of full-band count rate derived in bins of 
Alog(Count Rate) = 0.6 from stacking analyses, for all AGNs, all galaxies, and all sources (including both AGNs and galaxies), respectively. Horizontal dotted 
lines show the band ratios corresponding to given effective photon indexes, (b) Fraction of new sources as a function of full-band count rate for the main-catalog 
sources. The fractions are calculated in bins of Alog(Count Rate) = 0.6. [see the electronic edition of the Supplement for a color version of this figure.} 



the X-ray source positions in right ascension and declination 
and recorrelating with the optical sources. 

We present these 36 X-ray sources in Table 6 as a supple- 
mentary optically bright Chandra source catalog. The format 
of Table 6 is identical to that of Table 3 (see § 4.4 for the de- 
tails of each column). We note that the source-detection crite- 
rion is P < 0.1 for the sources in this supplementary catalog, 
as opposed to P < 0.004 for the main-catalog sources. Ad- 
ditionally, we set the multiwavelength identification-related 
columns (i.e.. Columns 18-22) to the WFI P-band matching 
results. 

5.2. Properties of Supplementary-Catalog Sources 

We show in Fig. 13(b) the positions of the 36 sources in 
the supplementary optically bright Chandra catalog, with the 
30 new sources shown as filled circles. These 36 supplemen- 
tary sources have /?-band AB magnitudes ranging from 15.1 
to 23.9. We show in Fig. 18(b) the P-band magnitude ver- 
sus the full-band flux for these 36 sources, with the sources 
being color-coded based on their likely types. For compari- 
son. Fig. 18(d) shows the IRAC 3.6 jim magnitude versus the 
full-band flux for these 36 sources. A total of 12 (33.3%) 



of these 36 sources are likely AGNs; 22 (61.1%) of these 
36 sources are likely galaxies and they all lie in the region 
expected for normal galaxies, starburst galaxies, and low- 
luminosity AGNs; 2 (5.6%) of these 36 sources are likely 
stars. The majority of these 36 supplementary sources ap- 
pear to be optically bright. X-ray faint non-AGNs (e.g., A03; 
Hornschemeier et al. 2003) as a result of our selection crite- 
ria, and thus they are not representative of the faintest X-ray 
sources as a whole. A total of 31 (86.1%) of these 36 sources 
have either spectroscopic or photometric redshifts. Of the 5 
sources that have no redshift estimate, 2 are bright stars with 
their redshifts set to -1 .000; the other 3 have their photometry 
severely affected by a nearby bright source, thus no redshift 
estimates were available. 

6. COMPLETENESS AND RELIABILITY ANALYSIS 

We performed simulations to assess the completeness and 
reliability of our main catalog; such practice has been com- 
mon among X-ray surveys (e.g., Cappelluti et al. 2007, 2009; 
Puccetti et al. 2009). 

6.1. Generation of Simulated Data 
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TABLE 6 

Supplementary Optically Bright Chandra Source Catalog 



X-ray Coordinates 



Detection Probability 



Counts 



No. 


0:2000 


^2000 


logP 


WAVDETECT 


Pos Err 


Off-Axis 


FB 


FB Upp Err 


FB Low Err 


SB 


SB Upp Err 


SB Low Err 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


1 ... 


03 31 44.64 


-27 45 19.4 


-1.9 


-5 


0.8 


10.10 


35.4 


19.2 


18.0 


25.8 


-1.0 


-1.0 


2... 


03 31 55.98 


-27 39 42.8 


-2.0 


-5 


1.1 


11.25 


24.2 


14.1 


12.9 


14.6 


8.1 


6.9 


3... 


03 31 56.42 


-27 44 11.4 


-1.9 


-5 


0.6 


8.19 


39.5 


21.5 


18.9 


16.6 


12.4 


9.8 


4... 


03 31 57.24 


-27 45 37.2 


-1.3 


-5 


0.8 


7.38 


64.8 


-1.0 


-1.0 


16.7 


11.5 


10.3 


5... 


03 32 07.63 


-27 49 27.2 


-2.3 


-8 


0.5 


4.63 


23.4 


11.2 


10.0 


12.2 


6.7 


5.5 



Note. — Units of right ascension are hours, minutes, and seconds, and units of dechnatlon are degrees, arcmlnutes, and arcseconds. Table 6 is presented in its entirety in the 
electronic edition. A portion is shown here for guidance regarding its form and content. The full table contains 79 columns of information for the 36 X-ray sources. 
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Fig. 17. — Average (i.e., stacked) band ratio in bins of (a) redshift (0 < 
Z < 1, 1 < z < 2, 2 < z < 3, and z > 3) and (b) X-ray luminosity [log(Lx) < 
41.5, 41.5 < log(Lx) < 42.5, 42.5 < log(Lx) < 43.5, 43.5 < log(Lx) < 44.0, 
and log(Lx) > 44.0] for new AGNs (filled circles), old AGNs (open circles), 
new galaxies (filled squares), and old galaxies (open squares). The median 
redshift or X-ray luminosity in each bin is used for plotting. The number of 
stacked sources in each redshift or luminosity bin is annotated. 
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First, we produced a mock catalog that covers the entire 
CDF-S and extends well below the detection limit of the 4 Ms 
mock 0.5-2 keV flux Hmits of (2-3) xlO"'^ 
Source coordinates were assigned using a 
recipe by Miyaji et al. (2007) to include reahstic source clus- 
tering. In this mock catalog, each simulated AGN was as- 
signed a soft-band flux that was drawn randomly from the 
soft-band log - log S relation in the AGN population syn- 
thesis model by Gilli, Comastri, & Hasinger (2007). Each 
simulated galaxy has a soft-band flux drawn randomly from 
the soft-band galaxy log - log S relation of the "peak-M" 
model by Ranalli, Comastri, & Setti (2005). The AGN and 



galaxy integrated fluxes match the cosmic X-ray background 
fluxes. The minimum soft-band fluxes simulated (« 3 x 10"'^ 
erg cm~^ s~^ for AGNs and « 2 x 10"'^ erg cm"^ s"^ for galax- 
ies) are a factor of « 3-5 below the detection limit of the cen- 
tral 4 Ms CDF-S (see § 7.2); inclusion of these undetectable 
sources simulates the spatially non-uniform background com- 
ponent due to undetected sources. The soft-band fluxes of the 
simulated AGNs and galaxies were converted into full-band 
fluxes assuming power-law spectra with F = 1 .4 and T = 2.0, 
respectively. 

Second, we constructed event lists from 54 simulated 
AClS-1 observations of the mock catalog, each configured to 
have the same aim point, roll angle, exposure time, and aspect 
solution file as one of the CDF-S observations (see Table 1). 
The MARX simulator was used to convert source fluxes to a 
Poisson stream of dithered photons, and to simulate their de- 
tection by ACIS. These event lists represent only events aris- 
ing from the mock point sources. 

Third, we extracted the corresponding background event 
files that are appropriate to the simulated source event files 
from the real 4 Ms CDF-S event files. For each real event 
file, we masked all the events relevant to the main-catalog and 
supplementary-catalog sources and then filled the masked re- 
gions with events that obey the local probability distribution 
of background events. The resulting background event files 
include the contribution (« 0.5%) of unresolved faint sources 
that was also present in the MARX- simulated source event 
files. To avoid counting the contribution of unresolved faint 
sources twice, we removed 0.5% of the events at random in 
each background event file and then combined it with the cor- 
responding source event file. Thus we produced a set of 54 
simulated ACIS-I observations that closely mirror the 54 real 
CDF-S observations. 

Finally, we obtained a simulated merged event file (i.e., 
sum of source and background events) following § 3.1, con- 
structed images from this simulated merged event file for the 
three standard bands following § 3.1, ran wavdetect on 
each simulated combined raw image at a false-positive prob- 
ability threshold of 10"^ to produce a candidate-list catalog 
following § 3.2, and utilized AE to perform photometry (and 
thus compute P values) for the sources in this candidate-hst 
catalog following § 3.2. 

6.2. Completeness and Reliability 

Our simulations allow us to assess the completeness and re- 
liability of our main catalog. Completeness is defined as the 
ratio between the number of detected sources (given a specific 
detection criterion P < Pq) and the number of input simulated 
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Fig. 18. — (Top) WFI R-band magnitude vs. full-band flux for sources in (a) the main catalog and (b) the supplementary optically bright catalog [note that the 
legend in Panel (b) appUes for all the panels in this figure]. Sources that are considered AGNs, galaxies, and stars are colored red. black, and blue, respectively. 
Open circles indicate AGNs/galaxies that were previously detected in (a) the LOS main catalog or (b) the LOS main or supplementary optically bright catalog; 
open stars in (a) indicate stars that were previously detected in the LOS main catalog; filled circles and stars indicate new AGNs/galaxies and stars, respectively. 
Arrows indicate limits. Diagonal lines indicate constant flux ratios between the WFI /?-band and the full band, with the shaded areas showing the approximate 
flux ratios for AGNs (dark gray) and galaxies (light gray). (Bottom) IRAC 3.6 fj.m magnitude vs. full-band flux for sources in (c) the main catalog and (d) the 
supplementary optically bright catalog. All the symbols are the same as those in Panels (a) and (b). The diagonal hnes indicate constant flux ratios between the 
IRAC 3.6 fim band and the full band. Note that several galaxies that have S-band detections were not detected in the IRAC 3.6 fj,m band, probably due to source 
blending in the IRAC 3.6 fim band and/or these galaxies being very blue systems, [see the electronic edition of the Supplement for a color version of this figure.} 



sources, above a specific source-count limit (this source-count 
limit applies to both the detected sources and the input sim- 
ulated sources). Reliability is defined as 1 minus the ratio 
between the number of spurious sources and the number of 
input simulated sources, above a specific source-count limit 
(again, this source-count limit applies to both the spurious 
sources and the input simulated sources). The top panels of 
Fig. 20 show the completeness and reliability as a function 
of the AE-computed binomial no-source probability P within 
the central 6 <6' area for the simulations in the full, soft, and 
hard bands, for sources with at least 15 counts and 8 counts. 
The bottom panels of Fig. 20 correspond to the case for the 
entire CDF-S field. 8 counts is close to our source-detection 
limit in the soft band. In each energy band, the completeness 
level for the case of 8 counts is, as expected, lower than that 
for the case of 15 counts, for both the central 6 <6' area and 
the entire CDF-S field; and the completeness level for the case 
of either 8 counts or 15 counts within the central 6 <6' area 
is higher than the corresponding completeness level in the en- 



tire CDF-S field. At the chosen main-catalog P threshold of 
0.004, the completeness levels within the central 6 < 6' area 
are 100.0% and 75.8% (full band), 100.0% and 94.1% (soft 
band), and 100.0% and 68.6% (hard band) for sources with at 
least 15 and 8 counts, respectively. The completeness levels 
for the entire CDF-S field are 82.4% and 49.3% (full band), 
95.9% and 63.5% (soft band), and 74.7% and 47.6% (hard 
band) for sources with at least 15 and 8 counts, respectively. 
The reliability level ranges from 99.2% to 99.8% for each en- 
ergy band and each source-count limit, which implies that, in 
the main catalog (i.e., the entire CDF-S field), there are about 
4, 4, and 3 spurious detections with > 15 counts in the full, 
soft, and hard bands, and about 4, 5, and 3 spurious detections 
with > 8 counts in the full, soft, and hard bands, respectively. 

We show in Fig. 21 the completeness as a function of flux 
under the main-catalog P < 0.004 criterion for the simulations 
in the full, soft, and hard bands. These curves of completeness 
versus flux derived from the simulations approximately track 
the normalized sky coverage curves (i.e., the curves of survey 
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Fig. 19. — Distributions of X-ray-to-optical flux ratio for (a) new AGNs 
(solid histogram) and old AGNs (dashed histogram) and (b) new galaxies 
(soUd histogram) and old galaxies (dashed histogram) with median flux ratios 
indicated by vertical lines. Shown in the insets are the distributions of R-band 
magnitude for new AGNs/galaxies (solid histograms) and old AGNs/galaxies 
(dashed histograms). Only sources with both full-band and i?-band detections 
are shown. 



TABLE 7 
Flux Limit and Completeness 



Completeness 


/o.5-8 keV 


/0.5-2 keV 


/2-8 keV 


(%) 


(erg cm~^ s~') 


(erg cm~^ s"') 


(erg cm~^ s~') 


90 


2.8 X l(r'5 


7.3 X IQ-'" 


4.0 X 10-'^ 


80 


1.3 X 


3.5 X ur'" 


2.0 X 


50 


3.0 X Kr'" 


7.8 X i(r" 


4.6 X 10"'* 


20 


1.1 X lO-'*" 


2.9 X 10-" 


1.8 X lO""" 



solid angle versus flux Umit; shown as solid curves in Fig. 21) 
derived from the real CDF-S data (see § 7.2). Table 7 gives 
the flux limits corresponding to four completeness levels in 
the full, soft, and hard bands, as shown as horizontal dotted 
lines in Fig. 21. 

7. BACKGROUND AND SENSITIVITY ANALYSIS 

7.1. Background Map Creation 

We created background maps for the three standard-band 
images as follows. We first masked the 740 main-catalog 
sources and the 36 supplementary catalog sources using circu- 
lar apertures with radii of 1.5 (2.0) times the w 99% PSF EEF 
radii for sources with full-band counts below (above) 10,000. 
Larger masking radii were used for the brightest sources 
(there are 3 main-catalog sources with full-band counts above 
10,000) to ensure their source photons were fully removed. 



Approximately 18.3% of the pixels were masked. By design, 
the background maps include minimal or no contributions 
from the sources in the main and supplementary catalogs; 
however, the background in the regions of a few extended 
sources (e.g., Bauer et al. 2002; LOS; A. Finoguenov et al. 
in preparation) will be slightly elevated. We then filled in the 
masked regions for each source with background counts that 
obey the local probability distribution of counts within an an- 
nulus with an inner radius being the aforementioned masking 
radius and an outer radius of 2.5 (3.0) times the « 99% PSF 
EEF radius for sources with full-band counts below (above) 
10,000. We summarize in Table 8 the background properties. 
We find our mean background count rates to be in agreement 
with those presented in L08. Our background is the sum of 
contributions from the unresolved cosmic background, parti- 
cle background, and instrumental background (e.g., Marke- 
vitch 2001; Markevitchet al. 2003). We do not distinguish 
between these different background contributions because we 
are here only interested in the total background. Even with a 
4 Ms exposure, the majority of the pixels have no background 
counts; i.e., in the full, soft, and hard bands, « 65%, 89%, and 
72% of the pixels are zero, respectively. 

7.2. Sensitivity Map Creation 

According to Table 4, the minimum detected source counts 
are w 11.4, 6.0, and 10.7 in the full, soft, and hard bands 
for the main-catalog sources, which correspond to full, soft, 
and hard-band fluxes of w 3.5 x 10"", 8.8 x lO^^^ and 6.4 x 
10"'^ erg cm"-^ s"', respectively, assuming that sources having 
these minimum counts are located at the average aim point 
and have a F = 1 .4 power law spectrum with Galactic absorp- 
tion. This calculation provides a measure of the ultimate sen- 
sitivity of the main catalog, which, however, is only relevant 
for a small central region near the average aim point. We cre- 
ated sensitivity maps in the three standard bands for the main 
catalog in order to determine the sensitivity as a function of 
position across the field. 

In the binomial no-source probability equation [i.e., eq. (1) 
in § 4. 1], we need to measure fisrc and 5ext to obtain the min- 
imum number of counts required for a detection (5), given 
the criterion of threshold = 0.004. We determined Bsrc in the 
background maps for the main catalog using circular apertures 
with 90% PSF EEF radii. Due to the PSF broadening with 
off-axis angle, the value of Bext has an off-axis angle depen- 
dency, i.e., the larger the off-axis angle, the larger the value 
of fiext- To follow the behavior of AE when extracting back- 
ground counts of the main-catalog sources, we derived the 
value of Bext as follows: for a given pixel in the background 
map, we computed its off-axis angle dp and set the value of 
fiext to the maximum Bext value of the main-catalog sources 
that are located in an annulus with the inner/outer radius be- 
ing 6p-0.25'/6p + 0.25' (note that the adopted maximum Bgxt 
value corresponds to the highest sensitivity). Given the com- 
puted Bsrc and B^xt, we numerically solved eq. (1) to obtain 
the minimum counts 5 (in the source-extraction region) re- 
quired for detections under the main-catalog source-detection 
criterion P < 0.004. We then created sensitivity maps for the 
main catalog using the exposure maps, assuming a F = 1.4 
power-law model with Galactic absorption. The above pro- 
cedure takes into account effects such as the PSF broadening 
with off-axis angle, the varying effective exposure (due to, 
e.g., vignetting and CCD gaps; see Fig. 2), and the varying 
background rate across the field. There are 1 1 main-catalog 
sources lying w 1-9% below the derived sensitivity hmits. 
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Fig. 20. — (Top) Case of 6 < 6': completeness (solid and dashed-dot curves; left y-axis) and reliability (long dashed and short dashed curves; right y-axis) 
as a function of Pq, the AE-computed binomial no-source probability threshold chosen for detection, for the simulations in the full, soft, and hard bands, for 
sources with at least 15 counts (red solid and long dashed curves) and at least 8 counts (blue dashed-dot and short dashed curves), respectively. Note that the 
short dashed curves overlap almost exactly along the long dashed curves in some cases (e.g., top-left and top-right panels). The vertical dotted lines indicate the 
chosen main-catalog source-detection threshold of Pq = 0.004. (Bottom) Same as top panels, but for the entire CDF-S field. \see the electronic edition of the 
Supplement for a color version of this figure.] 



TABLE 8 
Background Parameters 





Mean Background 


Total Background*^ 


Count Ratio'' 


Band (keV) 


(count pixel"' )" (count Ms"' pixel"' )" 


(10^ counts) 


(Background/Source) 


Full (0.5-8.0) 


0.482 0.252 


33.3 


10.6 


Soft (0.5-2.0) 


0.119 0.063 


8.2 


4.3 


Hard (2-8) 


0.363 0.178 


25.1 


20.5 



" The mean numbers of background counts per pixel measured from the background maps (see §7.1; note 
that a pixel has a size of 0.492"), which were not corrected for vignetting or exposure-time variations. 

The mean numbers of counts per pixel divided by the mean effective exposures (i.e., 1.909 Ms, 1.877 
Ms, and 2.040 Ms for the full band, soft band, and hard band, respectively) that are measured from the 
background maps (see § 7.1) and exposure maps (see § 3.1), respectively; these calculations take into 
account the effects of vignetting and exposure-time variations. 
' Total numbers of background counts in the background maps. 

** Ratio between the total number of background counts and the total number of detected source counts. 



i.e., 5 sources in the full band, 6 sources in the soft band, 
and none in the hard band, probably due to background fluc- 
tuations and/or their real T values deviating significantly from 
the assumed value. 

We show in Figure 22 the full-band sensitivity map for 
the main catalog. It is apparent that higher sensitivities are 
achieved at smaller off-axis angles. The sal arcmin^ re- 
gion at the average aim point has mean sensitivity limits of 
w 3.2 X 10"'^ 9.1 X 10-'^ and 5.5 x 10"'' erg cm'^ s"' for 
the full, soft, and hard bands, respectively. 



7.3. Sensitivities of and Prospects for Longer Chandra 
Exposures 



To investigate the improvement in sensitivity due to addi- 
tional exposure, we also created exposure maps, background 
maps, and sensitivity maps for the 1 Ms and 2 Ms CDF-S and 
simulated exposure maps, background maps, and sensitivity 
maps for the 8 Ms CDF-S. We followed the procedure de- 
tailed in § 3. 1 to create exposure maps for the 1 Ms and 2 Ms 
CDF-S. We simulated the 8 Ms CDF-S exposure maps by ro- 
tating the 4 Ms CDF-S exposure maps 90 degrees clockwise 
about the average aim point (this rotation approach simulates 
the variations of roll angles between observations) and then 
adding the rotated 4 Ms exposure maps to the real 4 Ms ex- 
posure maps. We followed § 7.1 to create background maps 
for the 1 Ms and 2 Ms CDF-S. To obtain the 8 Ms CDF-S 
background maps, we first simulated a set of 4 Ms CDF-S 
background maps by filling in each pixel in a simulated back- 



28 




Flux (erg cm ^ s ') 



Fig. 21. — Completeness as a function of flux under the main-catalog 
P < 0.004 criterion for the simulations in the full (blue filled circles), soft 
(green open diamonds), and hard (red open squares) bands, overlaid with the 
corresponding sky coverage curves (solid curves) calculated in § 7.3 and nor- 
malized to the maximum sky coverage (see the solid curves in Fig. 23). The 
dashed lines connect the corresponding cross points. The horizontal dotted 
lines indicate five completeness levels, [see the electronic edition of the Sup- 
plement for a color version of this figure.} 
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Fig. 22. — Full-band sensitivity map for the main catalog, created follow- 
ing § 7.2. The gray-scale levels, from black to light gray, represent areas 
with flux limits of < 4.0 X 10"", 4.0 X 10"" to 10"'^, 10"'^ to 3.3 X 10"'^, 
3.3 X 10"'* to 10"'^, and > 10"'^ erg cm"- s"', respectively. The regions 
and the plus sign are the same as those in Fig. 1 . 

ground map with background counts that obey the local prob- 
ability distribution of counts derived from the corresponding 
real 4 Ms CDF-S background map; we then rotated the simu- 
lated 4 Ms background maps 90 degrees clockwise about the 
average aim point and added the rotated 4 Ms background 
maps to the real 4 Ms background maps. We followed the 
procedure detailed in § 7.2 to create sensitivity maps for the 



1 Ms, 2 Ms, and 8 Ms CDF-S under the source-detection cri- 
terion P < 0.004, where we appropriately scaled the value of 
Bext that was derived when creating the 4 Ms sensitivity maps 
(i.e., scaling factors of 0.25, 0.50, and 2.0 were adopted for 
the 1 Ms, 2 Ms, and 8 Ms CDF-S, respectively). 

We show in Figure 23 plots of solid angle versus flux limit 
in the three standard bands for the 1-8 Ms CDF-S under 
the source-detection criterion P < 0.004. It is clear that, for 
each of the three standard bands, the quantitative increases in 
sensitivity are comparable between the cases of 1 to 2 Ms, 

2 to 4 Ms, and 4 to 8 Ms. To examine the improvement 
in sensitivity more clearly, we created sensitivity improve- 
ment maps by dividing the 1 Ms, 2 Ms, and 4 Ms sensitivity 
maps by the 2 Ms, 4 Ms, and 8 Ms sensitivity maps, respec- 
tively. We show in Figure 24 plots of solid angle versus mini- 
mum factor of improvement in sensitivity in the three standard 
bands between the 1 Ms, 2 Ms, and 4 Ms CDF-S and the 2 Ms, 
4 Ms, and 8 Ms CDF-S, respectively. Figure 24 only consid- 
ers the central 6* = 8' area, since such an area will be covered 
by any individual CDF-S observation. It is clearly shown in 
Fig. 24 that (1) for the three standard bands, the majority of 
the central CDF-S area generally has a factor of > V2 im- 
provement in sensitivity for each doubling of exposure time 
(note that V2= 1.414 corresponds to the background-limited 
case under the assumption, here inapplicable, of Gaussian 
statistics); and (2) among the three standard bands, the im- 
provement in sensitivity is most pronounced in the soft band 
for each doubling of exposure time, due to the fact that the 
soft band has the lowest background level (see, e.g.. Table 8). 
We note that, for each of the three standard bands during each 
doubling of exposure time, the improvement in sensitivity 
greater than a factor of 1 .5-1 .6 generally occurs in the ACIS -I 
CCD gap areas (see Fig. 2) where the improvement in expo- 
sure time is often greater than a factor of two. For the central 
w 100 arcmin^ area, the average improvement in sensitivity 
is typically a factor of 1.4-1.6 for each of the three standard 
bands, no matter which case of 1 to 2 Ms, 2 to 4 Ms, or 4 to 
8 Ms is considered. Based on the above analyses, we conclude 
that additional exposure over the CDF-S region, e.g., doubling 
the current 4 Ms exposure, will still yield higher sensitivi- 
ties in the central area of the field by a comparable amount 
to any previous doubling of exposure time (i.e., 1 to 2 Ms, or 
2 to 4 Ms). The faintest sources detected in an 8 Ms CDF-S 
should have full, soft, and hard-band fluxes of « 2.1 x 10"'^, 
6.0 X 10"'**, and 3.7 x 10"'^ erg cm~^ s"', respectively. Based 
upon the derived sensitivity maps and CXRB synthesis mod- 
els (e.g., Gilli, Comastri, & Hasinger 2007; Treister, Urry, & 
Virani 2009), a total of w 1000 sources, including « 120-130 
new AGNs and ss 90-100 new galaxies, are expected to be 
detected in an 8 Ms CDF-S. 

In addition to the improvements in sensitivity described 
above that would probe unexplored discovery space, signif- 
icant additional CDF-S exposure could greatly improve the 
X-ray spectra, light curves, and positions for the nearly 800 
known X-ray sources in our main and supplementary cata- 
logs. This would provide improved physical understanding of 
these sources; e.g., AGN content and luminosity, level and na- 
ture of AGN obscuration, shape of the X-ray continuum, and 
level of X-ray emission from X-ray binaries and supernova 
remnants. 



8. SUMMARY 
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Fig. 23. — Survey solid angle as a function of flux limit in the three stan- 
dard bands for the main catalog (shown as solid curves), determined follow- 
ing § 7.2. For comparison, the curves for the 1 Ms CDF-S (shown as dashed 
curves), the 2 Ms CDF-S (shown as dash-dot curves), and the simulated 8 Ms 
CDF-S (shown as dotted curves) are also plotted; these curves are calcu- 
lated/simulated consistently (see § 7.3). 
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Fig. 24. — Survey soUd angle within the central 6 = 8' area as a function of 
minimum factor of improvement in sensitivity in the three standard bands un- 
der the source-detection criterion P < 0.004. The improvement curves for the 
cases of 1 to 2 Ms, 2 to 4 Ms, and 4 to 8 Ms are shown as dash-dot, solid, and 
dotted curves, respectively; the curves for the first two cases are calculated 
using real data, while the curves for the third case are simulated (see § 7.3). 
The vertical dashed lines indicate a factor of \/2 improvement in sensitivity 
that is expected for each doubling of exposure time in a background-limited 
case under the assumption, here inapplicable, of Gaussian statistics. 

We have presented catalogs and basic analyses of X-ray 
sources detected in the deepest Chandra survey: the 4 Ms 
CDF-S. We summarize the most-important results as follows. 

1. The entire CDF-S consists of 54 individual observa- 
tions, with a summed exposure of 3.872 Ms and a total 
solid angle coverage of 464.5 arcmin^. 

2. The main Chandra source catalog contains 740 sources 
that were detected with wavdetect at a false-positive 
probability threshold of 10"^ and satisfy our binomial- 
probability source-selection criterion of f < 0.004; this 
approach is designed to maximize the number of re- 
liable sources detected. These 740 sources were de- 
tected in up to three standard X-ray bands: 0.5-8.0 keV 



(full band), 0.5-2.0 keV (soft band), and 2-8 keV (hard 
band). 716 (96.8%) of these 740 sources have multi- 
wavelength counterparts, with 673 (94.0% of 716) hav- 
ing either spectroscopic or photometric redshifts. 

3. The supplementary Chandra source catalog consists of 
36 sources that were detected with WAVDETECT at a 
false-positive probability threshold of 10"^ and satisfy 
the conditions of having 0.004 < P < 0.1 and having 
bright optical counterparts (R < 24.0). 

4. X-ray source positions for the main and supplemen- 
tary Chandra source catalogs have been determined us- 
ing centroid and matched-filter techniques. The ab- 
solute astrometry of the combined X-ray images and 
X-ray source positions has been estabUshed using a 
VLA 1 .4 GHz radio catalog. The median positional un- 
certainty at the w 68% confidence level is 0.42'70.72" 
for the main/supplementary Chandra source catalog. 

5. Basic analyses of the X-ray and optical properties of 
the sources indicate that they represent a variety of 
source types. More than 75% of the sources in the 
main Chandra catalog are likely AGNs. Near the cen- 
ter of the 4 Ms CDF-S (i.e., within an off-axis angle 
of 3'), the observed AGN and galaxy source densities 
have reached 9800!j]S! deg-^ and 6900!^^jJ° deg-^, re- 
spectively. The majority of the sources in the supple- 
mentary optically bright catalog are Ukely normal and 
starburst galaxies. 

6. A total of 300 main-catalog sources are new, compared 
to the 2 Ms main-catalog sources. Of the 300 new main- 
catalog sources, ~ 64% are hkely AGNs while « 35% 
are likely normal and starburst galaxies (the remain- 
ing pa 1% are likely stars), reflecting the rise of nor- 
mal and starburst galaxies at these very faint fluxes. In- 
deed, based on our source-classification scheme, galax- 
ies become the numerically dominant population of 
sources appearing at 0.5-8 keV fluxes less than w 10 
erg cm"-^ s"' or luminosities less than w 10'*'^ erg s"^ 
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Simulations show that our main catalog is highly reli- 
able (e.g., < 5 spurious detections are expected in the 
soft band) and is reasonably complete (e.g., the com- 
pleteness level for the soft band is > 94% for sources 
with > 8 counts in the central 6 <6' area). 

The mean background (corrected for vignetting and 
exposure-time variations) is 0.252, 0.063, and 0.178 
count Ms"^ pixel"' for the full, soft, and hard bands, 
respectively; the majority of the pixels have zero back- 
ground counts. 

The 4 Ms CDF-S reaches on-axis flux limits of 
w 3.2 X 10"'^ 9.1 X 10"'^ and 5.5 x lO"'^ ergcm'^ s"' 
for the full, soft, and hard bands, respectively, a factor 
of 1 .5-1 .6 improvement over the 2 Ms CDF-S. Another 
doubling of the CDF-S exposure time would still yield 
higher sensitivities in the central area of the field by a 
comparable amount to any previous doubling of expo- 
sure time, thus providing a significant number of new 
X-ray sources that probe the key unexplored discovery 
space. 
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The CDF-S source catalogs and data products provided 
by this paper will be beneficial to many ongoing and future 
studies; e.g., a search for a population of heavily obscured 
AGNs at intermediate redshifts (Luo et al. 2011), X-ray 
spectral constraints on heavily obscured and Compton-thick 
AGNs at high redshifts (Alexander et al. 2011; Gilli et al. 
20 11), derivation of X-ray number counts for different source 
types and the evolution of normal-galaxy luminosity functions 
(B. D. Lehmer et al., in preparation), and a study of extended 
sources (A. Finoguenov et al., in preparation). The CDF-S 
will continue to be a premiere deep- survey field over the com- 
ing decades; the CDF-S imaging and spectroscopic coverage 
are superb and continue to improve. For example, the Cosmic 
Assembly Near-IR Deep Extragalactic Legacy Survey (CAN- 
DELS)52 ^iii utilize 7/5rAVFC3 to image the GOODS-S; tiie 
highest sensitivity will be achieved with 5-orbit observations 
in the central region of this field. Together, deeper Chandra 
and multiwavelength data will be critical to allow comprehen- 
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